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ABSTRACT 


The  stress  fields  in  the  components  of  a  unidirectional  composite  due 
to  shrinkage  and  external  loads  are  computed  for  20  matrix/reinforcement 
combinations  having  various  volumetric  contents.  Further,  the  load  trans¬ 
missions  between  loaded  and  unloaded  fibers  are  formulated  as  three- 
dimensional  elasticity  solutions.  The  instability  problem  of  a  composite 
is  treated  by  both  the  static  and  the  energy  method,  resulting  in  critical 
loads  and  buckling  wavelengths  which  depend  on  material  constants  and 
geometries.  The  theoretical  results  are  in  good  agreement  with  experiments. 
The  work  reported  herein  encompasses  the  following  principal  areas: 

1.  Parametric  studies  (internal  stresses  and  displacements  com¬ 
puted  for  unidirectional  composites  composited  of  different 
combinations  of  matrices  and  reinforcements  and  different 
volumetric  contents) 

2.  Three-dimensional  load  transfer  among  loaded  and  unloaded 
fibers  in  a  matrix 

3.  Buckling  of  fibers  in  a  matrix  under  axial  load  as  an  elas¬ 
ticity  solution 

4.  Buckling  of  fibers  in  a  matrix  under  axial  load,  solved  with 
the  Ritz-Galerkin  method 

5.  Buckling  of  fibers  in  a  matrix  due  to  matrix  shrinkage 
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CHAPTER  i 


INTRODUCTION 


The  s trength-to-weight  efficiency  of  composite  materials  such  as  glass 
reinforced  plastics  has  been  demonstrated,  through  both  theoretical  and 
experimental  investigations,  to  be  superior  to  that  of  present-day  metallic 
materials.  While  the  application  of  these  materials  to  airframe  structures 
obviously  would  be  beneficial,  their  practical  utilization  has  been  limited 
because  reliable  structural  data  have  not  yet  been  developed,  and  because 
the  micromechanical  influence  of  the  composite  constituents,  along  with 
their  failure  initiation  and  crack  propagation,  has  not  been  fully  under¬ 
stood  . 

During  the  performance  of  Contract  DA  44-1 77-AMC -208 (T) ,*  mathematical 
relationships  were  derived  for  a  single  fiber  embedded  in  a  resin  cylinder 
for  the  case  of  static  loading  with  general  (mathematical)  boundary  con¬ 
ditions.  These  relationships  were  later  extended**  to  encompass  the  uni¬ 
directional  multifiber  composite  subjected  to  forces  of  thermal  contraction 
during  the  cure  cycle,  as  well'  as  to  external  loads. 

This  report  describes  a  continuation  of  these  efforts  to  define  the 
mechanical  behavior  of  fiber  reinforced  plastic  composites  in  order  that 
optimized  materials  and  structural  concepts  can  be  developed  for  airframe 
components  of  US  Army  aircraft.  Work  was  devoted  to  the  following  areas: 

1.  Internal  stresses  and  displacements  for  unidirectional  com¬ 
posites  composed  of  different  matrix/reinforcement  combinations 
having  different  volumetric  contents 

2.  Three-dimensional  load  transfer  among  loaded  and  unloaded 
fibers  in  a  matrix 

3.  Buckling  of  fibers  in  a  matrix  under  axial  load  as  an  elas¬ 
ticity  solution 

4.  Buckling  of  fibers  in  a  matrix  under  axial  load,  solved  with 
the  Ritz-Galerkin  method 

5.  Buckling  of  fibers  in  a  matrix  due  to  matrix  shrinkage 


*  See  USAAVLABS  Technical  Report  65-58. 

**  Contract  DA  44-177-AMC-320(T) ;  USAAVLABS  Technical  Report  66-62. 


CHAPTER  2 


SUMMARY  OF  RESULTS 


The  three-dimensional  analysis  of  a  multifiber  composite,  such  as  that 
depicted  in  Figure  1, 


Figure  1.  General  Arrangement  of  a 
Multifiber  Composite 

reveals  stresses  which  may  not  be  easily  anticipated  by  simple  reasoning. 

For  instance,  under  an  external  compression  load  in  the  direction  of  the 
fibers,  a  tensile  stress  is  generated  perpendicular  to  the  fiber  at  cp  =  0° 
and  becomes  compressive  at  cp  =  30°.  For  the  tension  case,  these  stresses 
are  reversed.  Therefore,  radial  tension  will  exist  in  the  bond  regardless 
of  the  direction  of  the  axial  load.  For  a  200,000-psi  axial  load,  for 
example,  these  stresses  are  between  2000  and  4000*  psi  and,  therefore,  a 
weakening  of  the  bond  between  fiber  and  resin  occurs. 

One  would  expect  the  shrinkage  stresses  in  the  resin  to  be  universally 
tensile;  in  the  radial  direction  of  the  interface,  however,  these  stresses 
are  compressive  where  the  fibers  are  closer  together  (cp  =  0°),  and  are  ten¬ 
sile  where  the  fibers  are  further  apart. 

One  would  expect  the  stress  in  the  fiber  direction  due  to  resin  shrink¬ 
age  to  be  compressive.  However,  in  high-density  fiber  composites,  part  of 

*  In  a  64%  glass  reinforced  epoxy  composite. 
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this  stress  in  the  fiber  is  tensile  due  to  Poisson's  effect  and  the  high 
radial  compressive  stress  where  the  fibers  were  closer  together.  This  stress 
seems  to  increase  stability  in  closely  packed  composites.  In  an  irregularly 
packed  composite,  it  would  tend  to  prebuckle  the  reinforcement. 

The  instability  solution  indicates  that  the  force  required  to  buckle  a 
fiber  in  a  multifiber  composite  is  much  higher  than  expected  based  on 
simpler  analyses.  One  would  be  tempted  to  conclude  that  a  densely  packed 
composite  fails  under  loads  other  than  buckling. 

The  buckling  wavelength  was  also  an  outcome  of  this  investigation.  It 
is  remarkable  that  the  wavelength  obtained  from  an  analysis  with  interfacial 
shear  taken  into  account  and  an  analysis  where  no  interfacial  shear  was  con¬ 
sidered  yield  very  similar  results.  From  this,  one  might  conclude  that  the 
buckling  wavelength  in  a  composite  is  independent  of  the  inter  facial  shear. 

By  comparing  the  two  instability  analyses,  however,  the  buckling  force  was 
found  to  be  very  dependent  on  the  interfacial  shear. 

In  general,  the  results  of  the  parametric  studies  and  the  instability 
analysis  complement  each  other,  leading  to  the  same  conclusions. 

On  the  basis  of  the  stresses  computed  during  this  program,  it  has 
been  determined  that  a  composite  loaded  in  compression  can  have  a  higher 
fiber  density  than  one  loaded  in  tension. 
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CHAPTER  3 


ANALYTICAL  WORK 


PART  I  -  PARAMETRIC  STUDIES 

This  study  was  made  in  order  to  investigate  the  stress,  strain,  and 
displacement  fields  in  the  fiber  and  matrix  for  composites  under  external 
and  residual  shrinkage  loads.  The  following  parameters  were  varied:  fiber 
and  matrix  modulus,  fiber  radius,  fiber  length,  and  fiber  volumetric  content. 
Combinations  of  the  above  parameters  were  picked  in  order  to  provide  data  on 
some  presently  used  fiber  and  resin  materials  in  composites.  From  the 
analytical  work,  a  computer  program  was  written  in  order  to  obtain  stress, 
strain,  and  displacement  values  in  the  composite.  A  discussion  of  the 
analytical  work  and  the  computer  program  is  given  in  Appendix  I. 

In  general,  the  objectives  were  to  predict  the  weaknesses  in  a  compos¬ 
ite  which  would  be  created  by  residual  shrinkage  loads  and  external  loading, 
and,  if  possible,  to  deduce  simplified  equations  for  the  pertinent  stress, 
strain,  or  displacement  fields.  Figures  2  through  28  demonstrate  that  these 
curves  could  be  fitted  by  fairly  simple  equations  which  are  functions  of, 
say,  cp  or  z  ,  but  to  include  the  stress  level  in  such  simplified  equations 
and  determine  how  it  varies  as  a  function  of  fiber  and  resin  modulus  and  of 
Poisson's  ratio- and  fiber  volumetric  content  proved  to  be  very  difficult. 

Table  I  shows  the  various  parametric  combinations  which  were  chosen  in 
the  study.  The  cases  in  Table  I  are  comparable  to  boron/epoxy  and  E-glass/ 
epoxy  composites  under  residual  shrinkage  and  external  loading.  The  elastic 
moduli  of  boron  and  E-glass  were  taken  as  60  x  10s psi  and  10  x  106psi, 
respectively.  The  epoxy  modulus  was  varied  for  the  E-glass  case  to  obtain 
more  pronounced  effects  on  the  field  values  than  would  be  obtained  for  the 
boron  epoxy  cases. 

The  length  of  the  boron  fibers  were  varied  to  see  what  effects  would 
result  from  the  long  and  short  boron  fibers  .  The  values  for  the  fiber 
volumetric  content  were  picked  to  include  typical  composite  values  (VI  = 

64%  and  70%)  and  high  fiber  density  packing  (V*  =  80%).  The  radius  of  the 
E-glass  fiber  was  varied  because  the  E-glass  fiber  is  made  with  various 
radii,  while  boron  fiber  has  a  typical  radius  of  2  mils.  Poisson's  ratios 
for  all  cases  were  0.2  for  the  fiber  and  0.35  for  the  epoxy. 

Stress,  strain,  and  displacement  fields  were  obtained  for  each  case. 

The  most  revealing  effects,  however,  occurred  in  the  stress  fields.  Since 
it  is  also  experimentally  possible  to  measure  the  stresses,  the  study  con¬ 
centrated  on  the  changes  in  the  stress  field  due  to  variation  of  the  param¬ 
eters.  From  the  computer  results,  the  maximum  and  minimum  values  of  the 
six  stresses  occurred  at  the  fiber-resin  interface.  The  stresses  shown  in 
the  figures  are  therefore  the  stress  at  the  interface,  r  =  a  .  Thus,  these 
stresses  provide  an  indication  of  possible  weaknesses  and  areas  of  failure 
in  a  composite. 
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Residual  Shrinkage  Stresses  at 
Interface  for  Combinations  of 
Fiber  and  Resin  Moduli  at  V1  = 
Radial  Stress  versus  cp 


Residual  Shrinkage  Stresses  at 
Interface  for  Combinations  of 
Fiber  and  Resin  Moduli  at  V^-  = 
Resin  Hoop  Stress  versus  cp 


Figure  6.  Residual  Shrinkage  Stresses  at 
Interface  for  Combinations  of 
Fiber  and  Resin  Moduli  at  =  0.64 
Fiber  Axial  Stress  versus  9 


Figure  7.  Residual  Shrinkage  Stresses  at 
Interface  for  Combinations  of 
Fiber  and  Resin  Moduli  at  V*  =  0.64 
Resin  Axial  Stress  versus  9 


Figure  8.  Residual  Shrinkage  Stresses  at  Interface 
When  E*  =  60  X  10s  psi  and  =  0.38  x 
10s psi  for  Different  Fiber  Volumetric 
Content.  Radial  Stress  versus  9 


Figure  9.  Residual  Shrinkage  Stresses  at  Interface 
When  E^  =  60  X  106psi  and  =  0.38  X 
10s psi  for  Different  Fiber  Volumetric 
Content.  Fiber  Hoop  Stress  versus  9 


Figure  10.  Residual  Shrinkage  Stresses  at  Interface 
When  ■  60  x  10^  psi  and  E^  *  0.38  x 
lCPpsi  for  Different  Fiber  Volumetric 
Content.  Shear  Stress  versus  cp 


Figure  11.  Residual  Shrinkage  Stresses  at  Interface 
When  =  60  x  lO^si  and  *  0.38  x 
106psi  for  Different  Fiber  Volumetric 
Content.  Resin  Hoop  Stress  versus  cp 


Residual  Shrinkage  Stresses  at  Interface 
When  =  60  x  lC^ppi  and  *  0.38  x 
106psi  for  Different  Fiber  Volumetric 
Content.  Fiber  Axial  Stress 
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Figure  13.  Residual  Shrinkage  Stresses  at  Interface  When 
El  =  10  x  10s  psi  and  =  0.38  X  1CP  psi  for 
Different  Fiber  Volumetric  Content.  Radial 
Stress  versus  <p 


Figure  14.  Residual  Shrinkage  Stresses  at  Interface  When 
E1  =  10  X  10s  psi  and  E11  =  0.38  X  K^psi  for 
Different  Fiber  Volumetric  Content.  Fiber 
Hoop  Stress  versus  cp 
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Figure  15.  Residual  Shrinkage  Stresses  at  Interface 
VJhen  =  10  x  10®  psi  and  E^-*  =  0.38  X 
lO^psi  for  Different  Fiber  Volumetric 
Content.  Shear  Stresses  versus  cp 


Figure  16.  Residual  Shrinkage  Stresses  at  Interface 
When  E^-  =  10  x  10s  psi  and  E*^  =  0.38  x 
10s  psi  for  Different  Fiber  Volumetric 
Content.  Resin  Hoop  Stress  versus  cp 
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Figure  19. 


Stresses  at  Interface  Due  to 
External  Load  for  Combinations 
of  Fiber  and  Resin  Moduli  at 
Vi  *  0.64.  Fiber  Axial  Stress 
versus  z /t 
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Figure  20. 


Stresses  at  Interface  Due  to 
External  Load  for  Combinations 
of  Fiber  and  Resin  Moduli  at 


0.64. 
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Resin  Axial  Stress 


Figure  21. 


Stresses  at  Interface  Due  to  External  Load 
When  E*  =  60  x  K^psi  and  E^  =  0.38  x  ltfpsi 
for  Different  Fiber  Volumetric  Content. 

Radial  Stress  versus  cp 


Figure  22. 


Stresses  at  Interface  Due  to  External  Load 
When  E1  =  60  x  10s  psi  and  E11  ■  0.38  x  lO^si 
for  Different  Fiber  Volumetric  Content. 

Fiber  Hoop  Stress  versus  cp 
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Figure 


Figure 


25.  Stresses  at  Interface  Due  to  External  Load  When 
E*  =  10  x  106psi  and  E**  =  0.38  x  106psi  for 
Different  Fiber  Volumetric  Content.  Radial 
Stress  versus  cp 


Different  Fiber  Volumetric  Content.  Hoop  Stress 
versus  cp 


Figure  27.  Stresses  at  Interface  Due  to  External  Load  When 
E*  =  10  x  10®psi  and  E-^  =  0.38  X  lO^psi  for 
Different  Fiber  Volumetric  Content.  Fiber  Shear 
Stress  versus  9 


Figure  28.  Stresses  at  Interface  Due  to  External  Load  When 
E*  =  10  x  10s  psi  and  E^  =  0.38  x  106psi  for 
Different  Fiber  Volumetric  Content.  Resin  Hoop 
Stress  versus  <p 
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TABLE  I 


CASES  USED  FOR  THE  PARAMETRIC  STUDY 


Par ameters 


General 
State  of 
Stress* 

Fiber 
Radius , 
10"3  in. 

Fiber 
Length , 
in . 

Epoxy 
Modulus , 
106  psi 

D 

2 

J 

0. 

38 

0 

.5 

10.0 

1 

j 

L 

, 

2 

5 

1. 

0 

5 

.0 

0. 

38 

I 

J 

2 

5 

1. 

0 

5 

0 

I 

1. 

0 

Fiber 
Modulus , 
10^  psi 


Fiber 
Vol . 
Content 
1 


60 


10 


64 

70 

80 

64 


70 

75 

80 

64 

70 

80 

64 


70 

80 

64 

64 


S  = 


resin  shrinkage  which  includes  both  cure  shrinkage  and 
differential  thermal  contraction. 


L  =  constant  load  applied  at  the  ends  on  the  resin. 


Also,  the  computer  results  indicate  that  the  stress  field  away  from  the 
ends  of  the  fiber  has  little  dependence  on  the  z  coordinate.  In  Figures  2 
through  7,  then,  the  stresses  which  are  not  plotted  as  functions  of  z/i, 
can  be  considered  as  constant  away  from  the  fiber  ends  with  respect  to  the 
z  coordinate. 

Residual  Shrinkage  Load 

From  the  cases  shown  in  Table  I,  the  effects  of  changing  the  following 
parameters  can  be  investigated  for  the  residual  shrinkage  load.  In  all 
these  cases,  the  fiber  shrinkage  was  zero  and  the  resin  shrinkage  was  l7o. 


1.  Changes  in  the  Combinations  of  Fibers  and  Resin  Modulus  at  =  0.64 

The  three  combinations  of  fiber  and  resin  moduli  are  shown  in  Table  II. 


TABLE  II 

COMBINATIONS  OF  FIBER  AND  RESIN  MODULI 


The  change  from  a  fiber  modulus  of  60  x  lO^psi  to  10  x  10s psi  did  not 
have  much  effect  on  the  shrinkage  stresses.  However,  changing  the  resin 

„  r-  I  II 

modulus  from  0.38  X  lOPpsi  to  1  x  lOP  psi  increased  the  stresses,  a-,1  ,  ct-,-,  , 
II  il  il 

and  a..  ,  by  a  factor  of  1  divided  by  0.38.  In  other  words,  when  the  ratio 

11  I 

of  resin  to  fiber  moduli  is  small,  the  change  of  residual  stresses,  a  , 

II  II  11 

a  ,  and  <j„„  ,  is  proportional  to  the  change  in  resin  modulus.  See  Fig- 

11  LL  I  II  I 

ures  2  and  3.  Three  other  stresses,  »  anc*  cr22  ’  are  s^own  *-n 

Figures  4  and  5. 

Figures  6  and  7  indicate  a  possible  ±20%  variation  in  the  axial  re¬ 
sidual  stress  in  both  the  resin  and  the  fiber. 

2.  Changes  in  the  Fiber  Volumetric  Content 

Figures  8  through  16  show  the  effects  on  the  stresses  due  to  fiber 
volumetric  content.  Figures  8  and  13  reveal  a  possible  weakness  in  the 
composite.  The  radial  stress  along  the  interface  has  positive  values.  If 
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this  stress  were  large  enough,  the  fiber-resin  interface  could  fail  and  the 
composite  effectiveness  would  be  lowered.  Also  the  shear  stress  at  the 
interface  (Figures  10  and  15)  would  additionally  stress  the  interface  bond. 
It  is  also  noted  that  increasing  the  fiber  volumetric  content  in  general 
increases  the  maximum  stresses  at  the  interface. 

The  fiber  axial  stress,  shown  in  Figure  12,  indicates  another  composite 
weakness.  Since  the  residual  stress  in  part  of  the  fiber  is  positive  for 
V*  =  0.80  ,  the  fiber  will  fail  at  a  lower  stress  level  when  a  tensile  load 
is  applied  to  it.  From  Figure  12,  the  residual  positive  stress  will  become 
even  higher  for  higher  volumetric  content. 

3.  Changes  in  Fiber  Radius  at  =  0.64 

The  fiber  radius  was  doubled,  from  0.0025  to  0.005;  no  effects  within 
the  accuracy  of  the  calculations  could  be  seen  in  the  stress  fields. 

4.  Changes  in  the  Fiber  Length  at  V-^  =  0.64 

The  fiber  length  was  taken  at  10  inches,  3  inches,  and  0.5  inch. 

There  were  no  effects  within  the  accuracy  of  the  calculations  on  the  stress 
field. 

Externally  Applied  Load 

From  the  cases  shown  in  Table  I,  the  effects  of  changing  the  following 
parameters  can  be  investigated  for  an  externally  applied  load.  In  Fig¬ 
ures  17  through  28,  the  stress  is  divided  by  the  average  external  load. 

1.  Changes  in  the  Combination  of  Fiber  and  Resin  Moduli 

The  three  combinations  of  fiber  and  resin  moduli  for  fixed  fiber 
volumetric  content,  V*  *  0.64  ,  are  shown  in  Table  II.  In  general,  the 
stresses  in  the  plane  perpendicular  to  the  fiber  axis  were  not  affected 
much  (see  Figures  17  through  20).  However,  the  radial  stresses  at  the 
interface  would  indicate  weaknesses  in  the  bond.  Since  the  radial  stress 
goes  from  minus  to  plus  for  cp  =  0°  to  cp  =  30°  ,  there  will  be  tension  in 
the  bond  regardless  of  whether  the  axial  load  is  tension  or  compression. 

For  a  compression  load,  a£  <  0  ,  the  results  in  Figure  17  would  tend  to 
cancel  the  residual  stresses  in  Figure  2,  but  the  radial  stresses  for  a  ten¬ 
sile  load  will  add  to  the  residual  stresses  and  weaken  the  bond  even  further. 

The  axial  stresses,  an<2  cr^  >  for  the  three  combinations  were 

independent  of  cp  and  within  ±5%  of  the  stresses  calculated  by  the  follow¬ 
ing  simple  equations. 


21 


a 


I 

33 


I  _ 

_ al 

E1  V1  +  E11  (l  -  V1 


II 

CT33 


F11  “ 

E  a, 


E1  V1  +  E 


Hli 


-V1) 


These  equations  for  the  axial  stress  due  to  external  load  are  valid 
for  higher  fiber  volumetric  content,  as  will  be  described  in  the  next  sec¬ 
tion  . 


In  addition  to  the  radial  stress,  the  shear  stress  at  the  interface 
(see  Figure  18),  when  superimposing  the  residual  shrinkage  and  external 
load,  will  also  increase  the  total  shear  for  tensile  load  and  reduce  the 
shear  for  compressive  load. 

The  axial  stress  for  fiber  and  resin  is  shown  in  Figures  19  and  20 
respectively. 

2.  Changes  in  Fiber  Volumetric  Content 

The  values  for  the  fiber  volumetric  content  were  normally  0.64,  0.70, 
and  0.80.  However,  an  unusual  change  in  the  fiber  hoop  stress  occurred  when 
fiber  content  was  changed  from  0.70  to  0.80  (see  Figures  22  and  27).  As  a 
result,  an  additional  value  of  fiber  content,  0.75,  was  taken  to  determine 
whether  the  fiber  hoop  stress  change  was  abrupt  or  gradual.  From  Figure  22, 
this  stress  distribution  occurred  between  a  fiber  content  of  0.70  and  0.80, 
and  can  be  expected  to  be  more  pronounced  beyond  0.80  fiber  content. 

The  radial  stresses  from  external  load  were  similar  to  the  previous 
results.  The  maximum  stress  is  positive  and  the  minimum  is  negative,  and 
they  both  increase  in  magnitude  as  the  fiber  volumetric  content  increases. 

The  axial  stresses  were  of  the  same  magnitude  as  calculated  from  the 
equations  in  the  previous  section. 

3.  Changes  in  the  Fiber  Radius  at  =  0.64 

There  was  no  effect  on  the  stress  fields  due  to  a  change  in  the  fiber 
radius  from  0.0025  inch  to  0.005  inch  for  externally  applied  load  at 
constant  volumetric  content. 


PART  II  -  THREE-DIMENSIONAL  LOAD  TRANSFER  AMONG  THE  FIBERS  IN  A  MATRIX 


The  fact  that  uneven  loading  occurs  in  a  composite  material,  especially 
near  the  boundaries  and  near  microfractures,  makes  it  necessary  to  under¬ 
stand  the  load  transfer  among  the  filaments.  Since  an  arbitrary  arrangement 
would  create  insurmountable  mathematical  difficulties,  the  mathematical  model 
chosen  was  one  of  periodically  repeating  characteristics.  The  unidirec¬ 
tional  fibers  were  assumed  to  be  in  a  hexagonal  array,  with  every  third 
fiber  axially  stressed  by  an  external  load.  The  rest  of  the  fibers  were 
assumed  to  be  stressed  by  load  transmission  through  the  matrix,  which  was 
the  main  subject  of  this  portion  of  the  investigation. 


A.  General  Arrangements  and  Assumptions 

The  composite  is  considered  to  be  free  from  residual  stresses,  and 
axially  finite  but  laterally  infinite.  Figure  29  depicts  part  of  the  cross 
section  of  the  composite,  with  the  shadowed  circles  representing  the  exter¬ 
nally  stressed  fibers. 


Figure  29.  Geometric  Arrangement  of  Fibers  in  a 
Unidirectional  Multifiber  Composite 


Displacements  and  stresses  of  a  typical  hexagonal  segment  in  the  com¬ 
posite  will  be  determined  by  assuming  that  the  segment  before  deformation 
changes  size  yet  remains  hexagonal  after  deformation.  However,  an  un¬ 
deformed  cross-sectional  plane  of  the  segment  may  not  remain  plane  during 
deformation.  Therefore,  the  segment  is  treated  by  the  classical  theory  of 
elasticity  as  a  three-dimensional  body. 


Since  displacements  and  stresses  in  fibers  and  matrix  are  hexagonally 
symmetric,  a  triangular  prism  of  the  hexagonal  segment  with  its  cross  sec¬ 
tion  A  0  C  D  and  coordinate  systems  as  shown  in  Figure  30  will  be  studied. 


Figure  30.  Geometry  of  a  Typical  Segment 
in  the  Composite 


The  natural  way  to  attack  the  problem  is  to  divide  the  prism  into  loaded 
fiber  (domain  I),  matrix  (domain  II),  and  unloaded  fiber  (domain  III),  at 
their  interfaces.  The  elasticity  solution  derived  from  Papkovitch  functions 
is  then  applied  to  each  of  these  domains  by  determining  arbitrary  constants 
with  appropriate  physical  conditions.  The  present  domains  are  now  defined 
as 

Domain  I  (Loaded  Fiber): 

Osr^a  ,  Oscp^r  >  -l  <  z  £  i 


Domain  II  (Matrix): 
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Domain  III  (Unloaded  Fiber) : 

[  2  2 

[a  sin  0  +  (b  -  a  cosG)  I  £  r  s' 

0  <.  cp  < 

--C  5  Z  < 

The  quantities  a  ,  b  ,  h  ,  l  ,  and  0  are  defined  in  Figure  30. 
B .  Formulation  of  Physical  Conditions 
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The  symmetry  of  radial,  circumferential,  and  axial  displacements 
g.  (i  -  1,2,3)  with  respect  to  zx-  and  xy-planes  dictates  that  in 
loaded  fiber,  matrix,  and  unloaded  fiber, 
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where  the  superscript  k  refers  to  the  corresponding  domains  I, 


(1) 

II,  or  III. 
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The  periodic  characteristic  of  displacements  in  loaded  fiber,  matrix, 
and  unloaded  fiber  leads  to 


s£(r,cp,z)  -  §^(r,9+5,z)  (2> 

where  k  =  I, II,  and  i  =  1,2,3  . 
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(3)  are  expressed 
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Figure  31.  Periodic  Characteristic  of 
the  Unloaded  Fiber 
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These  values  are  given  by 


Equation  (2)  states  that  the  displacement  vectors  of  matrix  and  exter¬ 
nally  stressed  fiber  possess  sixfold  symmetry.  Equations  (3)  are  the  con¬ 
dition  that  radial,  tangential,  and  axial  displacements  of  the  unloaded 
fiber  with  respect  to  its  axis  are  identical  for  every  multiple  of  2tt/3 
of  the  polar  angle  8  . 

If  continuity  in  displacements  and  appropriate  stresses  cy 

at  the  interface  exist,  then 


§^(a,<p,z)  =  ^I(a,cp,z) 


CTli(a,<p,z)  =  aii(a,cp,z) 

^(W2)  "  5™(r3'’’3'z)  (5) 

where  i  =  1,2,3  . 
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For  a  given  value  of  ©  in  the  range  0  £  ©  £  -j  ,  the  values  of  r3 

and  cp3  calculated  by  equation  (7)  are  the  coordinates  of  a  point  at  the 
interface  between  domains  II  and  III.  Equations  (6)  are  the  continuity  con¬ 
dition  on  normal  stress  in  the  p  direction  and  circumferential  and  axial 
shearing  stresses  at  the  interface  between  domains  II  and  III. 

In  order  that  the  fibers  remain  in  the  regular  hexagonal  arrangement 
during  deformation,  the  displacement  normal  to  the  hexagonal  boundary  must 
be  a  constant,  denoted  by 


§^(h  seci)r,<p,z)  cosi|/  +  seci|i,cp,z)  sinijf  =  6 

k  =  II, III 


(8) 
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where 


(9) 


Note  that  the  vanishing  of  the  circumferential  displacements  §2  at  the 
lines  OC  and  OB  (Figure  30)  will  be  automatically  satisfied  once  the 
symmetry  condition  (1)  is  fulfilled. 

Now  that  the  hexagonal  boundary  is  a  symmetrical  plane,  shearing 
stresses  vanish  at  the  boundary.  Thus,  in  the  directions  of  the  boundary 
line  CD  and  the  z  axis, 

al^(h  seci|(,cp,z)  -  a^Ch  sect)r,cp,z) 


-  2  cr^O1  secij/,cp,z)  cot  2ijr  =  0 


k  1c 

a13(h  sec\)r5cp,z)  +  seci|/,cp,z)  tanijj  =  0 
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k  =  II,  III 


If  the  axial  stress  applied  at  the  ends  of  the  central  fiber  as  shown 
in  Figure  30(a)  is  denoted  by  <jq  ,  the  appropriate  condition  is 
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and  if  the  ends  of  other  fibers  and  matrix  are  free  from  axial  stresses, 
then 
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It  can  be  seen  from  the  hexagonal  symmetry  that  shearing  forces  acting 
at  the  ends  of  a  hexagonal  segment  automatically  vanish  as  soon  as  the  con¬ 
dition  (2)  is  fulfilled.  If  the  ends  are  required  to  be  free  from  shearing 
stresses,  then 


CT3i(r  ,cp,±-t)  =  0 

k  =  I, II, III  i  =  1,2 


(13) 


The  formulation  of  physical  conditions  is  now  complete.  In  total,  we 
have  54  conditions  to  be  satisfied  in  this  analysis.  The  solution  derived 
from  Papkovitch-Neuber  functions  for  the  deformation  of  each  domain  must 
satisfy  appropriate  conditions  as  described  above. 


C .  General  Displacements  and  Stresses 

An  elasticity  solution  derived  from  Papkovitch  functions  was  given  in 
a  previous  report  (Reference  1).  If  nonperiodic  terms  are  taken  into 
account,  the  solution  to  Laplace's  equation  V^Pj  =  0  can  be  written  as 
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j  =  0,3 


where  P0  is  a  scalar  Papkovitch  function  and  P3  the  component  of  the 
Papkovitch  vector  in  the  axial  direction,  and  where  a  >  0  >  y  »  6  >  £  > 

Q  ,  T\  t  and  jj,  are  arbitrary  constants.  The  solution  Pj  can  be  expressed 
in  terms  of  Bessel  functions  of  first  and  second  kinds  of  r  ,  hyperbolic 
functions  of  z  ,  and  trigonometric  functions  of  cp  .  However,  solution 
(14)  has  the  advantage  that  the  argument  of  modified  Bessel  functions 
In  (jir)  and  Kn  (pr)  is  real  in  this  work. 

The  solution  to  the  following  equations, 

V2  |p^  coscp  -  s:*-ncp|  =  0 

V2  |Pj  sincp  +  P^  cosepj  =  0  (15) 

can  be  written  in  a  similar  form,  and  the  Papkovitch  functions  Pf  and  Pz 
in  radial  and  circumferential  directions  can  be  determined  from  these  two 
solutions . 

Displacement  components  are  related  to  Papkovitch  functions  by 
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Thus,  the  general,  displacement  function,  expressed  in  terms  of  cylindrical 
polar  coordinates,  can  be  found  and  simplified,  as  shown  in  Appendix  IX.* 

Stresses  of  an  elastic  body  are  related  to  strains  by  generalized 

Hooke's  law  as  follows:  J 
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*  Equations  (146)  through  (148). 
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The  stresses  corresponding  to  the  above  displacement  function  are  com¬ 
puted  from  equations  (16),  (17),  and  (18).  The  equations  for  stresses  are 
given  in  detail  in  Appendix  II.* 


D.  Determination  of  the  Arbitrary  Constants 

Arbitrary  constants,  which  enter  through  equation  (14)  into  solu¬ 
tions  (16)  and  (17),  are  to  be  determined  by  appropriate  conditions  of 
each  domain.  Solutions  (16)  and  (17)  have  been  expressed  by  the  solu¬ 
tions  in  Appendix  II.**  To  satisfy  the  periodic  condition  (2)  of 
matrix  and  loaded  fiber,  we  specify,  after  substituting  solutions  (14) 
and  (15)  into  (16), 


n  =  6,  12,  18,  24,  ... 


(19) 


Displacements  (16)  specified  by  (19)  have  been  applied  by  Haener  and 
Ashbaugh  (Reference  2)  to  a  unidirectional  multifiber  composite  under  re¬ 
sidual  shrinkage  and  axial  load.  The  choice  of  (19)  results  in  eight  of 
the  integration  constants  in  the  displacements  becoming  zero  for  the  domains 
I  and  II  respectively.***  In  order  that  the  displacements  and  stresses 
remain  bounded  at  the  center  of  the  fiber,  the  coefficients  8jnk  ■'•n  e9ua“ 
tion  (14)  must  be  zero  for  the  solutions  pertaining  to  the  fiber.  Satisfying 


*  Equations  (149)  through  (154). 

**  Equations  (146)  through  (154). 

***  Equation  (156). 
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the  loading  condition  (1),  we  choose  in  the  expression  for  CT33 
Appendix  II*)  the  pi  so  that  cospi-t  =  0  ;  therefore, 


(see 


I  kn 
V  =  2 1 


(20) 


k  =  1,3,5  ... 

In  the  stress  expression  CTo^  for  the  fiber,  all  tt  e  terms  involving  cosz 
will  automatically  be  zero  by  reason  of  the  above  choice  (20).  However, 
the  terms  involving  sin  z  will  not  be  zero  and  will  contribute  a  loading 
at  the  end  which  is  a  function  of  r  .  Since  this  is  compatible  with  the 
boundary  condition  (11),  the  coefficients  of  sine  in  the  stresses  must 
be  zero.**  Therefore,  additional  relationship  between  these  coefficients 
can  be  obtained.*** 

Additionally,  the  unloaded  condition  (12)  of  domain  II  leads  to 


II  _  krr 
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k  =  1,3,5  .  .  . 


(21) 


and  to  the  discarding  of  six  further  integration  constants  and  a  relation 
between  four  additional  integration  constants. t 

The  periodic  condition  that  displacements  of  the  unloaded  fiber  are 
identical  for  every  multiple  of  2tt/3  with  respect  to  0  has  been  de¬ 
scribed  by  equation  (3).  The  displacements  in  domain  III  do  not  have  any 
periodic  propertytt  in  the  polar  angle  cp  ,  the  fourteen  terms  contained  in 
the  P  functions  in  (16)  and  (17)  being  characterized  by  the  factor 
sin  nep  or  cos  nep  for  domain  III. 


*  Equation  (154) . 

**  Equation  (154). 

***  Equation  (155) . 
t  Equation  (157) . 

tt  As  far  as  a  hexagonal  element  as  shown  i 
displacements  of  the  six  unloaded  fibers 
ever,  an  elasticity  solution  is  not  vali 


Figure  30(a)  is  concerned, 
possess  sixfold  symmetry.  How- 
for  a  discontinuous  medium. 
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The  resulting  solution  is' the  special  case  corresponding  to  the  vanish¬ 
ing  of  n  in  the  expressions  P  and  ^ j  in  the  solutions  (16)  and  (17). 
The  vanishing  of  the  axial  stress  at  z  =  in  equation  (12)  requires  at 

the  end  of  domain  HI 


III 
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(22) 


k  =  1,3,5 


On  the  basis  of  equations  (22)  and  (12),  conclusions  similar  to  that 
made  for  the  loaded  part  can  be  made  in  this  case.* 

It  can  be  seen  that  displacements  and  stresses  are  composed  of  products 
of  functions  of  r  and  orthogonal  functions  of  cp  and  z  .  The  functions 
of  (p  will  no  longer  be  orthogonal,  however,  as  displacements  (16)  are 
introduced  into  conditions  (3)  and  (8)  and  stresses  (17)  into  conditions 
(6)  and  (10).  This  arises  from  the  following  parameters: 

1.  The  polar  radius  r  is  related  to  the  polar  angle  cp  at  the 
hexagonal  boundary. 

2.  The  polar  coordinates  r  and  cp  at  the  interface  between 
domains  II  and  III  are  related  to  each  other  by  condition  (7). 

3.  Coordinates  r  and  cp  in  domain  III  are  related  by  con¬ 
ditions  (4)  . 

Because  of  the  existence  of  certain  relations  between  r  and  cp  ,  it 
is  impossible  to  determine  the  arbitrary  constants  in  solutions  (16) 
and  (17)  by  conditions  (3),  (6),  (8),  and  (10)  as  they  stand.  Therefore, 
the  following  approximation  is  developed  such  that  these  conditions  will  be 
satisfied  with  any  desired  or  sufficient  accuracy. 

In  the  design  of  a  fiber  reinforced  composite  structure,  a  high  racio 
of  fiber  to  matrix  is  desired,  and,  hence,  the  angle  subtended  by  the  hex¬ 
agonal  boundary  of  the  matrix  toward  the  z-axis  is  small  in  a  cross-sectional 
plane.  For  example,  this  angle  is  about  5°  for  the  volume  of  fibers  in  a 
composite  being  657«.  Therefore,  the  hexagonal  boundary  of  the  matrix  (not 
fibers)  can  be  replaced  by  a  circular  cylindrical  surface  with  its  axis 
coinciding  with  the  z-axis  and  with  its  radius  being 

r  =  r  [ 2h  +  V4h2  +  (b  -  2a)2]  (23) 

m  4  L 


*  Refer  to  equation  (160). 
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With  such  a  substitution,  conditions  (10a)  and  (10b)  become 


au(V‘P>z]  '  0 

'’“(v?-2)  "  0  (24) 

The  approximation  to  conditions  (3),  (6),  (8),  and  (10)  will  now  be 
developed. 

Based  on  the  periodic  condition  (3)  of  the  unloaded  fiber,  average  dis¬ 
placement  components  of  a  volume  element  in  the  direction^  at  0=0  and 
0  =  rr/2  and  in  the  axial  direction  are  respectively  equal  to  those  of  a 
corresponding  element  in  the  directions  at  0  =  2tt/3  and  0  =  7tt/6  and 
in  the  axial  direction. 


t  1,2,3, . . .m^ 
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where  is  the  number  of  volume  elements  of  domain  III,  nd  and  Vj 

are  two  corresponding  equal  elements  as  shown  in  Figure  31  (.shadowed  areas 
with  the  same  axial  length). 

According  to  conditions  (6a)  through  (6c),  the  average  displacement  of 
the  matrix  at  an  element  of  the  intcrfacial  surface  between  domains  II  and 
III  is  equal  to  the  corresponding  average  displacement  of  the  unloaded 

fiber. 


i  =  1,2,3, ...m2 


which  are  the  corresponding  averages  of  the  displacement  components  along 
the  directions  of  x- ,  y- ,  and  z-axes  respectively,  and  in  which  is 

an  elemental  area  of  the  interfacial  surface  with  m2  being  the  number  of 
elemen  ts . 

From  conditions  (6d)  through  (6f),  two  averages  of  the  corresponding 
components  of  the  resultant  force  produced  by  normal  and  shearing  stresses 
at  an  element  of  the  interfacial  surface  between  domains  II  and  III  are 
equal  to  each  other. 
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rr 


A. 

1 


II 

rll 


r3 ’^3  5  z I  cos 


cp3  +  0jcoscp3  +  a 22  r3  >cPi  >z  |  s in | cp^  +  ej  sincp3 


a  12  ( r3  ,Cp3  ’ Z  )  sin(2  <P3  +  0) 


a  cos0)  d0  dz 


^“(w)  C0S(CP3  +  e)cos<P3  +  CT22I(r3’Cp3,Z)  sin(cP3  +  e)  sincp3  “ 


CT12I  (r3  ,CP3  ’  Z  j  sin  2  |^3  +  e) 


a  cos0  >  d0  dz 


all(r3,C(33’Z)  COS  ^3  +  0]sincP3  +  CT22  (r3  ,c^3  ’ Z)  sin(cf)3  +  0)  coscP3  + 


II 

f12 


( r3  ,Cp3  5  Z )  cos[2  ^3  +  0] 


a  s in0  >  d0  dz 


A. 

l 


all  r3’^3’Z]  C0S  ^3  +  0  sincP3  +  cr**1  r3’cP3>0]  Sin(cp3+  0  coscp3  + 


III  / 


CJ12I  (r3  ’^3  ’  Z)  cos(2  <P3  +  e) 


a  sin0  )  dQ  dz 


JJ 


A  *"• 

i 

rr 


|o  J3(r3,cp3,z)  cos |93  +  0)  -  CT23 (r3 ,Cp3  ’ ^  )  sirl(cp3  +  0]  ad0  dz 


JJ 


A. 

1 


Cr13I(r3’Cp3 jZ)  C0SK  +  9)  ‘  a23I(r3’CP3’Z)  sinK  +  6 


ad0  dz  (27) 


^>2j3>  ■  •  •  nin 


which  are  force  components  along  x-,  y-,  and  z-axes  respectively. 
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According  to  condition  (8),  the  average  of  the  displacements  normal  to 
the  hexagonal  boundary  must  be  equal  to  6  for  any  element  of  the  boundary 
plane 

CC  ~i 

1  f  k  k  2 

—  §^(h  secijf,cp,z)  cost  +  sect,cp,z)  sinf  h  sec  t  dijf  dz  =  6  (28) 


where  Bi  is  an  elemental  area  of  the  boundary  plane  of  the  matrix  or  the 
unloaded  fiber.  The  number  m3  of  the  elements  of  the  matrix  boundary  is 
not  necessarily  equal  to  that  of  the  fiber  boundary. 

Based  on  condition  (10),  the  shearing  forces  acting  at  an  elemental 
area  of  the  hexagonal  boundary  must  vanish. 


to 


^(h  sect ,cp,z)  -  (M  secij(,cp,z)  simjf  cost  - 


sect  ,cp,z)  (sin2?  “  cos2^J  h  sec  ♦  dz  *  0 


li 


sect  ,cp,z)  cost  +  o'23^ O’1  sect,cp,z)  sint 


h  sec  t  dt  dz  =  0  (30) 


i  =  1,2,3, . . .m, 


which  are  shearing  forces  along  the  boundary  line  C  D  (Figure  30)  and  the 
z-axis  respectively. 

Applying  a  similar  approximation  to  condition  (13),  the  vanishing  of 
the  shearing  force  which  acts  on  a  small  element,  R^,  of  the  end  section 
area  requires,  in  the  directions  of  x-  and  y-axes, 


Jr  ><t,<P’ 


,1)  cose p  -  cr.,-(r,cp,4)  sincp  rdcp  dr  =  0 


-O  sincp  +  Oo9(r,cp,i)  coscp  rdcp  dr  =  0 


k  -  I, II, III  i  =  1,2,3, . . .m. 
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It  is  seen  that  equations  (25)  through  (31)  can  be  reduced  to  condi¬ 
tions  (3),  (6),  (8),  (10),  and  (13)  when  the  regions  in  which  integrations 
in  these  equations  are  carried  over  become  infinitesimally  small.  There¬ 
fore,  these  conditions  are  replaced  by  the  corresponding  equations. 

If  two  functions  are  defined  as 


f  (z)  = 


z  (0  <  z  <  l) 

21  -  z  U  <  z  <  21) 


g(z) 


C z2  (0  <  z  <  O 

l-z2  +  hi  z  -  2t2  (1  <  z  <  21) 


(32) 


then  the  variables  z  and  z^  contained  in  the  P  functions  in  displace¬ 
ments  (16)  and  in  the  CT;Lj  of  stresses  (17)  can  be  represented  by 


1 

•  — r  sinjj,z 
k 


2 

z 


1 

— ^  COS|J,Z 

k 


(33) 


The  introduction  of  solutions  (16)  into  (25),  (26),  (27),  and  (28), 
and  (17)  into  (5d,e,f),  (27),  (24),  and  (31)  leads  to  a  system  of  linear 
algebraic  equations  in  terms  of  arbitrary  constants  and  the  unknown  6  . 

The  infinite  series  in  these  equations  are  to  be  truncated  in  numerical 
computation.  If  the  number  of  n  values  is  denoted  by  n  ,  the  number  of 
k  values  in  the  solutions  for  domains  I  and  II  by  Tq.2  ,  and  the  number 
of  k  values  in  the  solution  for  domain  III  by  k3  ,  then  the  total  number 
of  unknowns  is 


9  X  k12  X  n  +  6  x  1^2  +8x^  +  8  Xn  +  12 


The  total  number  of  equations  to  be  satisfied  is 


8  X  k.  .  X  n  +  5  X  k,  0  +  7  x  n  +  3m,  +  6m„  +  4m.  4-  6m.  +  6 

12  12  12  3  4 
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Thus,  infinite  series  existing  in  this  work  without  any  recurring  relation¬ 
ship  between  their  coefficients  are  to  be  truncated  or  the  numbers  ,  m2, 
m3,  and  m^  are  to  be  chosen  in  such  a  way  that  the  number  of  unknowns  is 
equal  to  the  number  of  equations.  However,  the  order  of  the  square  matrix 
is  extremely  large.  To  reduce  the  order,  the  matrix  is  subdivided  into  a 
number  of  rectangular  arrays,  and  each  array  in  turn  is  a  matrix.  After 
some  laborious  computation,  the  maximum  order  of  submatrices  yields 


Tc12  y  n  +  E  +  n  +  2 


or 


4  x  £3  +  4 

In  this  study,  symmetry  condition  (1)  of  a  hexagonal  segment,  periodic 
condition  (2)  of  matrix  and  loaded  fiber,  continuity  condition  (5)  at  the 
interface  between  matrix  and  loaded  fiber,  and  end  conditions  (11)  and  (12) 
have  been  exactly  satisfied,  while  the  other  physical  conditions  are  approx¬ 
imately  fulfilled.  It  is  expected  that  the  vanishing  of  shearing  stresses 
(10a)  and  (10b)  at  the  hexagonal  boundary  of  the  matrix  will  be  satisfied 
with  sufficient  accuracy  and  that  the  vanishing  of  shearing  stresses  (13a) 
through  (13d)  at  the  ends  of  matrix  and  loaded  fiber  and  the  constant  dis¬ 
placement  requirement  (8a)  at  the  hexagonal  boundary  of  the  matrix  will  be 
approximately  fulfilled  for  a  few  values  of  and  a  •  The  price  paid 

for  increasing  accuracy  in  the  fulfillment  of  other  conditions  is  that  the 
order  of  the  matrix  must  increase. 

The  variation  in  stresses  of  a  long  composite  with  every  other  fiber 
loaded  is  expected  to  be  small  in  the  middle  portion  of  a  hexagonal  segment, 
and,  hence,  the  corresponding  elements  chosen  for  computing  efficiency 
should  be  larger  than  those  close  to  the  ends.  The  transfer  of  loads  in 
the  present  case  possibly  is  an  end  effect;  that  is,  the  transmission  of 
the  largest  portion  of  loads  from  the  externally  stressed  fiber  in  a  hex¬ 
agonal  segment  to  six  unloaded  fibers  through  the  matrix  occurs  in  the 
neighborhood  of  the  ends. 


PARI  III  -  BUCKLING  OF  A  FIBER  IN  A  FINITE  ELASTIC  MATRIX 
UNDER  AXIAL  COMPRESSION 

This  study  was  concerned  with  the  determination  of  the  critical  load 
of  a  fiber  embedded  in  a  soft  elastic  matrix  subjected  to  axial  compression. 

The  finite  length  composite  cylinder  is  assumed  to  be  free  from  initial 
stresses.  When  the  load  increases  incrementally  from  zero,  the  fiber  remains 
straight  and  the  composite  is  under  compression  without  bending.  When  the 
load  continues  to  increase  and  reaches  a  certain  value,  the  fiber  deflects 
laterally.  At  this  value,  the  compressive  force  is  called  the  Euler  criti¬ 
cal  load.  Later  on,  the  deflection  increases  rapidly  with  a  small  increase 
of  the  applied  compression.  Eventually  the  fiber  buckles  in  a  wavy  shape 
and  loses  its  natural  function,  the  transmission  of  compressive  forces. 
Therefore,  the  critical  compressive  load  is  very  important  in  the  design 
of  such  composite  structures. 

An  analytical  method  has  been  developed  by  Sadowsky  and  Hussain  (Refer¬ 
ence  3)  to  determine  the  thermal  critical  load  of  an  infinite  fiber  bonded 
to  an  infinite  domain  of  matrix  without  mechanical  loading.  The  matrix  has 
been  treated  by  the  linear  theory  of  elasticity  as  a  three-dimensional  body, 
and  the  fiber  by  one-dimensional  elasticity.  The  method  of  approach  is 
reasonable.  In  the  present  work,  it  additionally  accounted  for  the  moment 
produced  by  axial  shear  at  the  interface  between  fibers  and  matrix  such  that 
this  shear  influences  the  critical  load.  It  should  be  pointed  out  that  the 
total  axial  shearing  force  at  the  interface  vanishes,  but  not  the  produced 
momen  t  s . 

The  method  of  approach  to  the  present  problem  is  similar  to  that  devel¬ 
oped  by  the  above  authors.  Instead  of  using  Boussinesque-Papkovitch  poten¬ 
tial  functions,  the  equilibrium  equations  in  cylindrical  coordinates  for  the 
deformation  of  matrix  are  directly  solved  in  this  work.  Equations  of  equi¬ 
librium  for  the  critical  load  of  the  fiber  are  in  a  simple  manner  derived 
by  the  one-dimensional  nonlinear  theory  of  elastic  stability,  based  on 
statics  consideration  different  from  that  in  the  report  by  Sadowsky  and 
Hussain.  The  contribution  of  axial  shear  at  the  interface  to  the  equation 
of  moment  equilibrium  is  taken  into  account.  The  critical  compressive  load 
of  the  composite  cylinder  corresponding  to  the  buckling  of  the  fiber  found 
in  this  study  is  a  function  of  elastic  constants,  fiber  radius,  and  outer 
radius  of  the  composite. 


General  Description 


Consider  a  composite  cylinder  of  finite  length  L  ,  fiber  radius  a  , 
and  outer  radius  b  of  the  composite  with  two  coordinate  systems  as  shown 
in  Figure  32 . 
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Figure  32.  Geometry  of  a  Composite  Cylinder 


Experiment  shows  that  the  buckling  mode  of  a  fiber  in  a  low-modulus 
matrix  subjected  to  axially  compressive  load  is  approximately  a  plane  sine 
curve.  If  stresses  of  the  matrix  at  the  buckling  of  the  fiber  are  within 
its  elastic  range,  it  will  be  appropriate  to  treat  the  matrix  by  the  linear 
theory  of  elasticity  as  a  three-dimensional  body  and  the  fiber  by  the  non¬ 
linear  theory  of  elastic  stability  as  a  slender  bar  subjected  to  compression 
and  appropriate  lateral  loading.  If  the  buckling  wave  is  chosen  in  the  yz- 
plane  (Figure  33)  and  if  the  deflection  of  the  fiber  in  the  y-direction  is 
denoted  by  v  ,  then  the  buckling  wave  of  the  fiber  can  be  written  as 


v  =  v  srn#z 
o 


where  vq  is  the  amplitude  of  the  sine  curve  and  where 


being  the  number  of  half  wavelength  along  the  z-axis. 
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Figure  33.  Fiber  at  the  Buckling 


Displacements  and  stresses  given  by  (16)  and  (17)  for  load  transfer 
through  fibers  in  a  matrix  (Part  II  of  this  chapter:  Three-Dimensional 
Load  Transfer  Among  the  Fibers  in  a  Matrix)  are  suitable  to  the  deformation 
of  the  matrix  under  consideration  by  appropriate  choice  of  the  Cartesian 
coordinate  system.  However,  arbitrary  constants  in  these  expressions  are 
difficult  to  determine  in  the  present  case.  To  simplify  the  problem,  the 
approach  described  in  Part  I  of  this  chapter  (Parametric  Studies)  is 
applied . 

Based  on  the  three-dimensional  theory  of  elasticity  for  a  close- 
packed,  fiber  composite  under  compression,  stresses  in  the  plane  per¬ 
pendicular  to  fiber  axis  were  vanishingly  small  compared  to  the  average 
axial  stress  for  the  epoxy-fiberglass  composite. 

It  may  be  assumed  in  the  present  analysis,  without  introducing  appre¬ 
ciable  error  in  the  results, that  the  stresses  developed  during  buckling  at 
the  interface  between  the  fiber  and  matrix  are  produced  by  bending  only. 

The  displacements  and  stresses  of  the  matrix  are  given  in  Appendix  III 
to  this  report. 

Equations  of  Equilibrium  of  the  Fiber  and  Its  Critical  Compressive  Load 

Consider  an  element,  AB  ,  of  the  fiber  with  length  dz  before  bending. 
It  is  suppositioned  that  this  element  has  already  been  under  compression. 
After  bending,  AB  will  displace  to  A'B1  . 

I; 

. 

When  some  lateral  deflection  is  produced  in  the  compressed  bar,  there 
is  some  change  in  compression,  but  a  more  detailed  investigation  by  S.  Timo¬ 
shenko  (Reference  4)  shows  that  this  change  is  negligible.  Since  the 
compressive  load  is  not  changed,  the  fiber  axis  will  be  under  the  same  axial 
force  as  that  before  bending, and  there  will  be  no  change  in  the  length  along 
the  neutral  plane.  During  bending,  axial  forces  N  ,  transverse  shearing 
forces  Q  ,  bending  moments  M  at  the  ends  of  this  element,  lateral 
forces  S  per  unit  length,  and  axial  shearing  forces  T  per  unit  length 
are  shown  in  Figure  34. 
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Figure  34.  Geometry  of  the  Fiber 


The  equation  of  equilibrium  on  the  moment  in  the  x-direction  requires 


f;  +  Q  +  Mi  -  0  (36) 

where  the  last  term  is  the  contribution  of  moment  per  unit  length  pro¬ 

duced  by  axial  shear  at  the  interface  between  fiber  and  matrix. 

If  the  radius  of  curvature  of  the  element  dz  is  denoted  by  R  ,  then 
the  result  of  two  transverse  shearing  forces  is  in  the  negative  direction 
of  the  z-axis,  with  the  magnitude  being  Qdz/R,  and  the  result  of  two  axial 
forces  is  Ndz/R  in  the  negative  direction  of  the  y-axis.  Therefore,  the 
equilibrium  equations  on  the  force  components  along  y-  and  z-axes  are  given 
by 


^  . 

dz 


0 


(37) 


dN 

dz 


+ 


T 


a 

R 


0 


(38) 


where  the  nonlinearity  has  been  introduced  by  considering  the  influence  of 
the  deflection  v  and  where  the  angle  of  rotation  of  the  element  has  been 
assumed  to  be  small. 


The  moment  M  produced  by  axial  shear  <7i3(a,cp,z)  in  the  x-direc- 
tion,* 

2  I  CT13  (p  2  ,CfJ>  Z)  sincp  ^ 
o 


M  =  2a‘ 

s 


2  o 

TT  a  v  g  cosaz 
0  Id 


(39) 


and  the  component  of  the  moment  produced  by  g^  at  the  interface  in  the 
y  direction  automatically  vanishes. 

The  lateral  force  S  at  the  interface  yields,  by  using  (196), 


'  2tt 


S  = 


all(P2,Cp,Z)  sincP  +  CTi2(P2,Cp,Z)  COSCp 


ad(p 


o  ,  o  , 

tt  a  v0  <7^  +  <j  I  sinctz 


(40) 


and  the  force  component  produced  by  stresses  g^  and  g-^  at  the  inter¬ 
face  in  the  x-dire'ction  automatically  vanishes. 

The  axial  shearing  force  T  at  the  interface  is  from  (196) 


■2tt 


T  = 


13 


|p2>cp,zj  adcp  =  0 


(41) 


Using  (34),  (39),  (40),  and  (41)  and  the  relations 


1 

a-* 

R 


d2v 

dz2 


(42) 


M  = 


E1  I 


g  4  I  d  v 
-  4  a  -E  — 

dz 


(43) 


*  See  equation  (196). 
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from  which  cv  or  the  ratio  of  number  of  half  wavelength  to  fiber  length 
can  be  computed.  The  compressive  force  corresponding  to  this  value  of 

oi  is  the  critical  load  of  the  fiber  denoted  by  P^r  .  The  corresponding 
compressive  force  applied  to  the  ends  of  the  matrix  is  determined  from  the 
condition  that  the  axial  displacement  of  the  matrix  is  the  same  as  that  of 


47 


the  microfiber  at  the  ends  of  the  composite.  If  this  load  is  denoted  by 
PH,  we  have  from  the  above  condition, 


The  total  compressive  load,  P*  ,  applied  at  the  ends  of  the  composite  is 
then  given  by 

P*  =  P1  +  P11  (55) 

cr 

which  corresponds  to  the  buckling  of  the  fiber. 

The  constant  axial  displacement  at  the  ends  of  the  composite  is  not 
necessary  to  require  the  vanishing  of  the  right-hand  side  of  (193c).  The 
axial  displacement  produced  by  bending  is  of  influence  with  axial  stresses, 
but  not  the  axial  force  acting  on  a  cross  section  because  of  geometrical 
symmetry  of  the  fiber  with  respect  to  the  neutral  plane.  Actually,  the 
axial  displacement  produced  by  bending  at  the  beginning  of  the  fiber  buck¬ 
ling  can  be  neglected,  since  the  length  of  the  fiber  axis  does  not  change 
by  bending  and  the  angle  of  rotation  is  small.  If  this  is  done,  the  cor¬ 
responding  change  in  other  equations  is  to  set  the  right-hand  side  of  (195) 
to  zero. 

It  is  seen  from  (52)  that  (i=l,2,3)  is  the  contribution,  to  the 

critical  load,  of  interfacial  pressure  and  circumferential  and  axial  shearing 
stresses  at  the  interface  between  fiber  and  matrix.  The  critical  load  P1 
and,  hence,  P^  and  P*  are  simply  functions  of  elastic  constants  of  fiber 
and  matrix,  fiber  radius,  and  outer  radius  of  the  matrix.  Thermal  buckling 
load  of  an  inf ini te- length  fiber  bonded  to  an  infinite  domain  of  the  matrix 
without  mechanical  loading  and  critical  load  of  a  f ini te- length  fiber  in 
lateral  infinite  matrix  under  axial  compression  can  be  treated  as  special 
cases  of  the  present  result  by  performing  the  integrals  in  (177),  (181), 
and  (185)  analytically. 
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PART  IV  -  ENERGY  SOLUTIONS  OF  THE  INSTABILITY  PROBLEM 


The  problem  of  the  buckling  of  microfiber  embedded  in  a  supporting 
matrix,  under  axial  compression  and  under  internal  shrinkage  load,  was 
analytically  investigated  with  energy  methods  in  five  different  cases. 

The  models  include  a  single-fiber  and  multifiber  reinforced  composite  of 
finite  dimension,  and  a  single  finer  in  an  infinite  matrix.  The  multifibers 
under  internal  shrinkage  loads  are  treated  with  a  simplified  assumption, 
allowing  the  use  of  the  Lagrange  mu  1  *"i " 1  ’  er  method  of  vari,.  '  luii. 

The  Ritz-Galerkin  method  minimizing  total  potential  energy  was  used  in 
deriving  the  equation  for  obtaining  the  critical  load  of  the  fiber.  In  the 
present  analysis,  total  potential  energy  contains  the  strain  energy  due  to 
the  bending  and  extension  of  the  fiber,  plus  the  energy  of  the  interfacial 
pressure  and  longitudinal  shear  applied  by  the  matrix  to  the  fiber  less  the 
work  done  by  external  force.  The  stress  field  of  the  binder  was  determined 
by  solving  a  plane  elasticity  problem.  However,  it  was  solved  two- 
dimensionally  as  well  as  three-dimensionally  for  the  case  of  the  multifiber 
reinforced  composite  of  finite  size.  The  fiber  and  matrix  elasticity  are 
well  matched  at  the  interface  by  imposing  boundary  conditions  requiring  the 
continuity  of  stresses  or  displacements.  Finally,  the  smallest  buckling 
criterion  was  found  by  variation  with  respect  to  the  wave  parameter. 

To  formulate  the  problem,  the  following  assumptions  were  made: 

1.  Both  materials  are  linearly  elastic,  isotropic,  and  homogeneous. 

2.  Elongation,  shear,  and  rotation  of  constituents  are  small  in 
comparison  with  unity. 

3.  The  deflection  curve  of  the  fiber  and  then  the  distribution  of 
the  interfacial  force  at  the  buckling  state  are  sinusoidal  in 
the  axial  direction. 

4.  Binder  modulus  is  much  smaller  than  fiber  modulus. 

5.  The  cross-sectional  area  of  the  fiber  remains  plane  and 
circular  after  loading. 

6.  During  buckling,  there  is  no  relative  slip  between  fiber  and 
matrix . 


49 


A. 


Buckling  of  Multifibers  in  a  Finite  Matrix 


1.  Buckling  Load  of  the  Fiber  (Ritz  Method) 


Figure  35.  Half  Critical  Wavelength 
Shown  on  One  Fiber 

Geometric  boundary  conditions  are,  at 

Z  =  0  ,  =  0 


and  at 

z  -  L  ,  =  0 


(56) 


(57) 


Taking  into  consideration  the  geometrical  boundary  conditions  of  the  fiber, 
we  can  assume  its  deflection  curve, 


§ 


I 

1 


(58) 


In  the  present  problem,  the  total  potential  energy  includes  strain 
energy  due  to  the  bending  and  extension  of  the  fiber,  plus  the  strain  energy 
due  to  interfacial  pressure  and  axial  shear  applied  by  matrix  on  the  fiber, 
less  the  work  done  by  external  force. 
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The  strain  energy  due  to  bending  of  the  fiber  is  shown  below  (see 
Figure  36) . 


Figure  36.  Bending  of  a  Fiber  Element 

The  strain  energy  due  to  the  extension  of  the  fiber  is  derived  as  follows 
(see  Figure  37) . 


Figure  37.  Schematic  Diagram  of  a  Fiber 
Element  During  Loading 


The  deformed  length  of  an  element  of  the  fiber  is 


ds  ~-\f[ dz  +  d  C3)2  +  jd  Cj)2 


=  dz 
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d  §; 
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1 
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l  dz  j 

Then  the  increase  of  strain  energy  due  to  the  shortening  of  the  fiber  is 


P  (ds  -  dz) 


rL 

r .  1 

,  I, 

2" 

[" 

fd^l 

dz  2 
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'  0 

- 

dz 


(60) 


The  strain  energy  due  to  the  interfacial  pressure  contributed  by  the  matrix 
is 


1  , 

2  Pln  *1  dZ 


1 

2 


kk,  dz 


(61) 


The  strain  energy  due  to  axial  shear  contributed  by  matrix  is  depicted  in 
Figure  38. 
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a 


h - L  - 

Figure  38.  Geometry  of  Two  Neighboring  Fibers 
in  a  Multifiber  Composite 

In  the  matrix,  shear  strain  is  expressed  by 


e 


II 

13 


dz 


Since  shear  strain  should  be  continuous  at  the  interface,  th 


B§“(a)  agJCa) 

-M  I  ■  =  ■!■■■  ■  I 

3z  Bz 

Also 

,  II  -II,  >,  II,.  ,,11/  s 
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3r  =  d  **  ~ d 


But 


then 


II 

c13 


(64) 


II 

£13 


Thus,  the  strain  energy  due  to  longitudinal  shear  of  the  matrix  on  the  fiber 
is 


The  external  work  by  prescribed  loading  during  buckling  is 

f  L 

w  -  PL  sj(0)  -  PR  sjd)  -  E3  p3  d* 

Jo 


-  PL  53(0>  -  PR  ?3a>  + 


Figure  39.  Schematic  Diagram  Showing  One 
Fiber  Under  External  Load 
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Substituting  equation  (58)  into  equation  (69)  and  making  use  of  integral 
identities,  we  get 


Minimization  of  total  potential  energy  with  respect  to  the  amplitude  of  the 
assumed  curve  gives  the  stationary  condition;  i.e., 


BT 


=  0 


ij 


j  =  1,2,3, 


From  equation  (70) 
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1  eI  ^  iL 


+ 


3 

TT 

2L 


b 


2-a2) 
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(71) 


I 

This  equation  either  yields  |lj  =  0  ,  which  gives  the  undesirable  case  of 
zero  deflection,  or  the  following  results  are  obtained: 


±E*  I1 


.2  3 

J  TT 


I 

cr 


0 
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which  gives  the  critical  load  for  the  fiber;  i.e., 


2  2 
n  tt 


,1  -I 


cr 


I  + 


2  2 
n  n 


k  +  n(b2-a2jG11  (l  +f) 


(72) 


where  k  has  the  dimension  of  pound  per  square  inch,  and  is  equal  to  the 
ratio  of  interfacial  force  per  unit  length  to  unit  lateral  deflection. 

Foundation  constant  k  for  the  multifiber  reinforced  matrix  is  derived 
as  shown  in  Figure  38. 

We  can  imagine  the  matrix  as  a  two-dimensional  plate  and  under  merely  the 
normal  interfacial  force  per  unit  length  p^(z)  .  The  expression  p-,(z) 
must  have  similar  distribution  as 
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’1°°  ■  I 
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Therefore,  the  boundary  conditions  are,  at  r  =  0  , 


and  at  r  =  2d  , 
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The  stress  components  are 


II 

ct33 
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Applying  boundary  conditions  (74)  to  (77),  we  have 
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Substitution  of  (85)  and  (86)  into  (87)  yields 


k  § 


II 

In 


rra 


co 


sh(^]  +  Bn 


cosh  1 2  raid 

L 


sinh 


1 2  rd) 

+  D 

I2  "nd-)  sinh 

2  nnd  \ 

l  L  / 

n 

L  / 

k  § 


II 

In 


na 


(89) 


Substituting  (85)  and  (86)  into  (88),  we  get 
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The  denominator  of  Bn  and  Dn  ,  solved  from  equations  (89)  and  (90),  U 


2  nnd  ]  ^ 


I  -  sinh 


,'„k2  f  2  -nnd 


(91) 


The  numerator  of  coefficient  Bn  ,  solved  from  equations  (89)  and  (90),  is 
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The  numerator  for  coefficient  Dn  ,  found  from  (89)  and  (90),  is 


II 

—  sinh  (  2_OTTd\ 
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Then  the  results  of  coefficients  are  as  follows: 
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(94) 
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By  integrating  equation  (100),  we  get,  after  extensive  analysis, 
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where 
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2.  Smallest  Buckling  Load  and  Buckling  Wavelength  of  the  Fiber 


From  equation  (72),  we  have 


Pi  .  xi  +  ^  k  +  V  lit  4 

L  n  tt  ‘ 


cr 


(103) 


Let  ~  =  a  .  Combination  of  equations  (102)  and  (103)  yields  the 


expression  for  critical  load 
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Since  the  buckling  load  (104)  depends  on  the  wavelength 


(105) 


(104)  must  be  minimized  with  respect  to  arn  to  obtain  the  critical  wave 
length  corresponding  to  the  smallest  critical  buckling  load  of  the  fiber 
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The  result:  of  (106)  is 
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F 1 

From  equation  (107),  we  can  determine  c*n  as  a  function  of  ■=—  , 

E11 

a  and  b  .  In  other  words,  we  obtain  the  critical  wavelength  L  -  — 

6  cr  cy  j 

as  a  function  of  the  material  constants,  the  diameter  of  the  fibers,  and 
the  volume  percentage  of  the  fibers.  Substituting  this  value  of  o-i  into 

equation  (104),  we  then  get  the  minimized  critical  buckling  load,  P1 

cr  min 
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3.  Alternative  Way  of  Finding  the  Strain  Energy  Due  to  Longitudinal 
Shear 
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1 
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II  II 
a13  e13 


dv 


Jv 


(108) 


The  expression  from  equation  (84)  has  been  introduced  into  (108). 

After  performing  the  indicated  operations,  the  results  of  strain  energy  due 
to  longitudinal  shear  have  been  used  to  replace  the  last  term  of  equation 
(103).  In  this  manner, then,  the  smallest  buckling  load  and  corresponding 
wavelength  were  obtained. 
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(109) 


Y  is  an  analytical  expression  given  in  Appendix  IV. A  (198a),  The  equation 
for  an  is  given  in  Appendix  IV. A  (200),  while  k  (a  kind  of  a  foundation 
constant)  is  given  in  equation  (102). 


B.  Buckling  of  a  Single  Fiber  in  a  Finite  Matrix 

1.  Constant  k  for  Single  Fiber  in  a  Finite  Matrix 


a 

1 
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b  1  0 

1  d  ■  b-a 

r  \ 

Figure  40.  Single  Fiber  Dimensions 


Assume  load  per  unit  length  at  the  interface, 


(110) 
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The  boundary  conditions  are 


at  r  =  0  , 


at  r  *  c  =  b-a  , 
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By  applying  the  solutions  for  stresses,  equations  (83)  and  (84)  the 
integration  constants  are  established  [Appendix  IV. B,  equations  (201) 
through  (204)].  These  constants  are  used  in  the  displacements  and  stresses 
to  express  the  strain  energy  during  loading 


1 

2 


dv 


(112) 


written  out  in  Appendix  IV. B,  equation  (214),  which  contains  the  strain 
energy  due  to  shear  at  the  interface.  The  critical  load  was  then  computed: 
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and 


The  expression  cvn  ,  X  in  equation  (113)  can  be  found  in  equation  (220) 
(215)  respectively,  while  dn  ■  and  . 
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C.  Buckling  of  a  Finite  Single  Fiber  in  an  Infinite  Matrix 


If  the  matrix  is  infinite,  then  b  =  »  .  Thus,  it  is  more  convenient 
that  the  solution  of  the  differential  equation  (79)  takes  the  following 

form 
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By  routine  procedures,  the  stresses  and  strains  are  found  from  the  potential 

$  .  * 

Pi 

k  =  4-  (115) 


where  p^  is  the  lateral  unit  interfacial  force  and  g.  is  the  lateral 
displacement.  The  minimum  critical  buckling  load  and  the  corresponding 
wavelength  are  derived  in  (236)  through  (237),  and  are  restated  here: 
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The  wavelength  corresponding  to  the  critical  load  is 
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The  numerical  results  of  the  buckling  wavelength  for  a  single  fiber 
coincide  with  test  results  obtained  from  experiment  performed  during  this 
contract.  In  Figures  42  through  60,  the  critical  wavelength  and  the 

E1 

buckling  loads  of  a  multifiber  are  plotted  as  a  function  of  _  ,  the 

E11 

fiber  radius  a  ,  with  volume  percentage  content  of  the  fiber  as  a  param¬ 
eter.  In  Figures  61  and  62,  the  wavelength  and  buckling  load  of  a  single 
fiber  are  plotted  versus  the  radius  of  the  fiber. 


*  Equations  (224)  through  (226)  and  (229). 
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Figure  49.  Critical  Buckling  Wavelength  of  the  Fiber 
versus  Radius  of  the  Fiber  (E  =  5 . 32  X 
107psi,  E11  =  3.8  X  105psi,  or 
El/El1  =  140) 
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Figure  51.  Critical  Buckling  Load  of  the  Fiber 
versus  the  Ratio  of  Young's  Moduli 
(a  =  1  X  10‘3in.) 


Figure  52.  Critical  Buckling  Load  of  the  Fiber 
versus  the  Ratio  of  Young's  Moduli 
(a  =  2  x  10_3in.) 
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Figure  57.  Critical  Buckling  Load  of  the  Fiber  versus 
Radius  of  the  Fiber  (El  =  2.28  X  lO^psi, 
EH  =  3,8  x  105psi,  or  EI/E11  =  60) 


Figure  58.  Critical  Buckling  Load  of  the  Fiber  versus 
Radius  of  the  Fiber  (E*  =  3.8  X  lO^psi, 

E11  =  3.8  X  105psi ,  or  EI/E11  =  100) 
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D .  Buckling  of  a  Multifiber  Composite.  Due  to  Resin  Shrinkage 

In  this  treatment,  the  energy  was  derived  for  a  fiber  resin  cell.  The 
equilibrium  equation  formulated  as  an  energy  variational  problem  derived  in 

Appendix  IV. D  is 
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E^A11  ge  \dz 
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(118) 


with  the  restraining  condition 


(1  +  e)  coscp  -  (1  +  g)>  dz  =  0 


(119) 


In  equation  (118),  the  first  brackets  contain  the  compression  energy 
in  the  axial  direction;  the  second,  the  bending  energy  of  the  resin  and  the 
fiber.  The  third  expression  in  (118)  is  the  shrinkage  energy  of  the  resin. 
The  right-hand  expression  in  (118)  is  the  energy  introduced  by  the  external 
force.  The  restraining  condition  (119)  is  derived  from  a  simplifying  assump¬ 
tion  that  the  longitudinal  change  of  the  matrix  is  equal  to  that  of  the 
fiber.  It  must  be  emphasized  that  in  this  manner,  the  shear  produced  in 
longitudinal  direction  during  buckling  is  neglected  to  simplify  the  resulting 
Lagrange  equation. 

The  Lagrange  equations  corresponding  to  (118)  are 


Pg  +  R  4-  \  coscp  =  0 

(120) 

Qcp"  +  \  (1  +  s)sincp  =  0 

(121) 
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The  factor  \  is  the  Lagrange  multiplier  expressed  in  Appendix  IV. D, equa¬ 
tion  (274).  From  equations  (120)  and  (121),  the  following  differential  for 
cp  is  obtained  : 


Qcp"  +  X  1 
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(122) 


The  characteristic  equation  of  (122)  and  the  derivations  necessary  are  given 
in  Appendix  IV. D*  and  are  restated  here: 


3 

9 
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(123) 


To  obtain  the  displacement  of  the  fiber  axis,  cp  has  to  be  substituted  into 
equation  (351)  of  Appendix  IV. D.  The  wavelength  is  obtained  by  equation 
(302) ,  or  by 
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(124) 


For  small  amplitudes  of  the  maximum  slope  cp  of  the  buckling  curve,  the 
Lagrange  multiplier  becomes,  for  the  first  approximation,** 

X  -  -  (R  +  0  P)  (125) 


This  is  a  reasonably  good  solution  for  the  eigenvalue  \  .  By  calculating 

the  second  approximation,  we  must  use  equation  (350)  in  Appendix  IV. D.  In 

this  case,  \  is  the  force  exerted  upon  the  fiber  by  the  resin. 

In  spite  of  the  fact  that  the  interfacial  shear  was  neglected,  the  re¬ 

sults  of  this  analysis  coincide  with  the  results  obtained  by  Rosen  (Refer¬ 
ence  6).  For  E1  =  10  x  lO^psi  and  E11  =  3.8  x  105psi  ,  where  radius 
(a)  =  2.5  X  10-3  ,  half  distance  of  the  fibers  (b)  =  2.976  x  10“3  ,  and  the 
shrinkage  0  =  1%. 


*  Equations  (275)  through  (388). 

**  Reference  equation  (393). 


The  values  in  (124)  and  (125)  are 


P  = 

R  = 

\  - 

Q  = 

The  computed  wavelength 


,I„I  .II^II 

A  E  +  A  E 

=  199.56 

„II  .11  „II 

E  A  g 

4.56  x  10 

-  (r  +  e  p)  - 

1.95 

E1!1  +  E11!11 

=  3.17  x 

z  =  0.114  in. 

w 


E.  Buckling  of  Multifibers  in  a  Matrix  Under  Axial  Load  With  the  Matrix 
Treated  as  a  Three-Dimensional  Cylinder 

The  critical  buckling  load  of  the  fiber  was  found  by  minimizing  total 
potential  energy  with  respect  to  the  amplitude  of  the  assumed  deflection 
curve  of  the  fiber'*  and  is  rewritten  as  follows: 

2  _I  TI  .  k  .  1,2  2\  II  /  a)2 

anE  1  +  —  +  "lb  "a  )G  l1  +  di 

a 


cr 


where  k  is  the  unit  interfacial  force  per  unit  lateral  deflection;  i.e., 


k 


is  inversely  proportional  to  the  critical  wavelength. 


Specifically , 


2tt 


O' 


j 


*  Equations  (70)  through  (72). 


82 


The  displacement  and  stress  fields  in  the  matrix  are  determined  by 
solving  a  three-dimensional  elasticity  problem  (see  Appendix  IV. E).  To 
obtain  critical  wavelength,  we  then  take  differentiation  of  P^r  with 
respect  to  cyn  (equation  437)  after  substitution  of  k  (equation  436) 
into  equation  (103),  The  expression  for  finding  the  value  of  a j  is 
given  in  equation  (438).  Therefore,  the  smallest  critical  buckling  load 
can  be  found  from  equation  (437)  with  the  introduction  of  critical  wave¬ 
length  . 
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PART  V  -  SHORT-FIBER  STUDIES 


During  the  contractual  period,  exploratory  work  was  performed  to  estab¬ 
lish  any  potential  existing  in  the  utilization  of  short  fibers  as  reinforcing 
material  in  composites.  Limited  experiments  were  made  to  determine  the 
elastic  moduli  and  strength  in  compression,  and  the  fatigue  behavior  of 
unidirectional  composite  in  tension.  The  specimens  were  comprised  of  con¬ 
tinuous  glass  fibers  and  short  boron  fibers  oriented  unidirectionally . 

The  results  of  these  exploratory  experiments  are  reported  in  Appendix  V. 

The  results  of  static  tests  in  compression  were  compared  with  results 
available  from  tests  made  on  pure  glass  fiber  and  pure  boron  fiber  composite, 
and  it  has  been  concluded  that  the  modulus  of  a  composite  can  be  adjusted  to 
any  level  between  that  of  the  basic  reinforcing  constituents  merely  by  vary¬ 
ing  their  volumetric  relationship  at  a  constant  resin  content. 

The  primary  purpose  of  the  fatigue  study  was  to  determine  whether  the 
effect  of  the  end  of  short  boron  fibers  is  one  of  destroying  the  matrix 
material  and  disintegrating  the  composite.  According  to  the  test  results, 
it  appears  that  this  might  indeed  be  the  case.  When  comparing  short  fiber 
boron  composite  with  continuous  glass  fiber  composite  and  with  7075  aluminum, 
the  very  limited  test  results  indicate  that,  although  the  short  boron  fiber 
is  somewhat  better  in  fatigue  behavior  than  aluminum,  it  is  inferior  to 
continuous  glass.  However,  an  additional  factor  which  must  be  considered  is 
that  due  to  the  higher  modulus,  the  deflections  under  load  are  smaller  than 
those  of  a  pure  glass  fiber  composite. 

An  analysis  of  elasticity  has  not  yet  been  applied  to  the  short  fiber 
composite.  However,  the  experimental  results  appear  to  justify  a  strict 
analytical  investigation  of  this  type  of  material. 
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APPENDIX  I 


PARAMETRIC  STUDIES 


ANALYTICAL  DISCUSSION 

The  following  conditions  and  assumptions*  pertain  to  the  computer 
program  which  was  developed  to  determine  the  three-dimensional  stress, 
strain,  and  displacement  fields  in  a  composite.  Also,  a  discussion  and 
listing  of  the  computer  program  are  given  at  the  end  of  this  Appendix. 

It  was  assumed  that  the  composite  consists  of  identical  hexagonal 
elements.  Each  element  has  a  centrally  located  fiber  surrounded  by 
resin  as  shown  in  Figure  62.  Because  of  the  symmetry  involved,  the 
following  domains  of  the  polar  coordinates,  r,  cp,  z,  are  sufficient  to 
determine  the  displacement  field  in  the  hexagonal  element. 


Figure  63.  Typical  Hexagonal  Element  of  a  Composite 


*These  conditions  and  assumptions  are  also  presented  in  Reference  1. 
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For  the  fiber: 


06  A.  «  cl 

0$  <p  *  3o° 

For  the  resin: 


(126) 


a.  *  a,  ^  cc^<p 
o  *  <P  £  30 

°*}s  * 

The  distance,  b,  is 


(127) 


given  in  terms 


of  the  fiber  volumetric  content,  V  : 


4- 


(128) 


At  the  interface,  r  =  a ,  the  displacements,  and  appropriate 

j  i —  ~  rtrtnhi'niinilfi  ' 


stresses,  au>  were  assumed  to  be  continuous 


Sift'S*  =  &*” 


L  -  IjZj3 


(129) 


At  the  hexagon  boundary  during  shrinkage  and  axial  loading  the  pi 
forming  the  hexagon  column  containing  Bfi  and  parallel  to  the  z  axis  or 
fiS  axis  remains  a  plane  because  it  must  be  shared  by  the  adjacent 

hexagonal  element. 

Relating  the  hexagonal  elements  in  Figure  62  to  a  rectangular  coordi¬ 
nate  svstem  x  y  and  z,  there  is  complete  symmetry  about  the  x  -  b 
plane!  Therefore,  the  shear  stresses  must  vanish  in  this  plane  so  that 


<r*  =  O 

-7 


°*r° 
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In  addition  to  the  vanishing  shear,  the  normal  displacement  of  BD 
relative  to  coordinates  at  point  A  must  be 


For  the  shrinkage  case,  e  / b  would  represent  the  total  lateral  com¬ 
posite  shrinkage.  For  the  loaSed  case,  it  would  represent  the  total  lateral 
expansion  or  contraction.  Specifying  the  normal  displacement  and  shear 
stresses  at  the  hexagon  boundary,  however,  has  introduced  another  unknown: 
e  .  An  additional  equation  is  obtained  by  considering  that  the  total  normal 
stress  on  the  hexagon  boundary  plane  x  =  b  in  Figure  62  has  to  balance 
the  external  side  loading.  Because  these  side  toads  are  zero  for  the  pres¬ 
ent  case,  the  normal  stress  must  satisfy  the  equation 

\  \  ^  ^  dl  -° 

The  conditions  for  a  ,  O  ,  and  §  on  the  hexagon  boundary  must  be 
satisfied  for  all  values  o¥  y  and  z,xsuch  that 

°  $  y  $ 

The  corresponding  values  in  cylindrical  coordinates  become 

A  =  ,  O  s  (pz  \  ,  OZ  ^  Sr  JL 

The  relation  between  the  stresses,  cr  ,  a  ,  and  a  ,  and  displacement, 
5  ,  with  reference  to  rectangular  Cartesian  coo?dimtesXand  the  stresses  and 
displacements  in  polar  coordinates  are 


<rKX  =  €Tn  -  <r/JL  OITL  (f 

a-  =  ^2  ( <r,.  -  <%z  )  z  (p  +  <r,z  go- a,  z  cp 

</ 

°T3  <P  -  <^3  S**'  <P 

^  ^  ■CSX1*'  <p  —  /a <p  (130) 
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The  boundary  conditions  at  the  ends  z  =  ±1  of  the  composite  are 

±1)  =  Of1 

<r3“  (a, $,±1)  =  cr/ 

Additionally,  the  following  shear  conditions  at  the  composite  ends  are 

fulfilled. 


5$ 

n 

ss 


tJi)  A  cLa,  ctcp  -  O 
(  A-j  <Pj  -  &)  A.  tXl/L  dcf)  -  O 

tJL)  /LclA-otty  -  o 
(Aj  <P,  ±^)  -^ok.dcf>  =  0 


For  the  shrinkage  case,  and  are  zero,  and  the  shrinkage 

coefficients  and  31  are  not  zero.  For  the  enci  loading  case,  3*  and 

3-LX  are  zero  and  oy-  and  have  constant  values.  By  superposition  of 

the  two  cases,  the  total  stress  is  obtained. 

The  materials  were  assumed  to  be  individually  isotropic  linear  materi¬ 
al.  The  general  assumptions  for  the  classical  theory  of  elasticity  were 
used.  From  the  many  possibilities,  the  approach  toward  determining  the 
solutions  of  the  Navier  equations  by  potentials  connected  with  the  names  of 
Bousslnesq,  Papkovich,  and  Neuber  has  been  selected.  These  general  and 
complete  solutions  by  potential  functions  for  three-dimensional  problems 
are  given  below  in  vector  form  to  facilitate  the  resulting  transformation  to 
cylindrical  coordinates. 
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and 


Si  *  <£2  *  ^3  S3  ]  =  tiu//  veetoe 

7^  =o  (134) 


Since  £.  and  &  are  functions  of  cp,  it  is  more  convenient  to  relate 

e,  and  e„  to  the  unit  vectors  e  and  e  associated  with  rectangular 

~1  , ,  ~2  ~x  ~y 

coordinates. 


e  (  =  gj, 

=  "^>  CXY&Cp 

$3  3  Sj  (135) 


Then,  by  substituting  (135)  into  the  first  equation  of  (134), 

V  [(Tfcycroff  -/? +  (/f  ryc*i*(f>  +  Pz  cmp')Sa  * 

7^  g  J  -  tiu.ll  uec 

d 

Since  e  ,  e  .  and  e  are  independent  of  position,  then 
'vx  ~z 

c^aro-cP  —  7^  <p  J  —O 

V2(  Tf  /&urt(p  -f  75  ocro-  —  0  (136) 

The  solution  of  equations  (136)  and  the  last  of  (134)  can  be  written 
in  terms  of  the  solution  to  Laplace's  equation. 


=0  (*'0,/A  3  (137) 

where 

4.0  5 

Lt  =  P  ocrxcp-  7J  f 

L  £  -  /yi^x.(p  +  7J  con-  cp 

43  =  ^  (138) 
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and  the  second  and  third  equations  of  (138)  can  be  solved  for  Pj_  and  P_ 
so  that  (i  =  0,  1,  2,  3)  can  be  written  in  terms  of  (i  =  0,  1,  2,  3). 


~P(  -  Lt  c-o-^cp  *  L-, 

~  L  i  OCsyl  (p  -h  l_ 


OC-vt,  Cp 

2 


3 


(139) 


The  strains  become 


and  the  stresses  are 


) 


cr 

4j 


v*  ET  r 


(140) 


(141) 


Therefore,  from  the  solution  of  (137),  the  Papkovitch  functions  can  be 
determined  by  (139),  then  the  general  displacements  by  (133),  and  then  the 
strains  and  stresses  by  (140)  and  (141)  respectively. 


DISPLACEMENT  AND  STRESS  EQUATIONS 

After  deriving  the  general  displacements  and  stresses,  equations  from 
(133)  through  (141),  some  of  the  symmetry  and  boundary  conditions  were  used 
to  reduce  the  equations  to  the  following  form.* 


*The  argument  of  all  the  Bessel  functions  is  a^r. 
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The  displacements  for  the  fiber 


5,X=  fl  21  ChL  ( * 


r  s+-4)/:T'n.  j-  _  «<*■  A.  7—  \  . 
'•4 nit  \  9(l-V~)  •*-*.+!  9(t-^r)  J  r 


C/W  C 4  -<ntT*)I°)  +^^hjTf]c^o(/e^ 


-  *  t-Zv^  K 

”>,a  *<>-»-) 


t-t.  £.U 


r  t-  +  j>r+n.  j-  £.r  +-fri.n  ^  ^ 

v/<  +<!-»*)  J”n-/  n£i'n/<  *{!•»-)  ■% h 


°C  Hr^r)n-In-  J 


/CW-W.  "tyf 


e1  -  y~  y  f/f  1 

S3  ~  £—  C—  L ^nrk  9(1-^“)  m'‘  '*■' 


^  1  0(^/1  — -  ^ 

'Zrnk  9 {i-'ST)  J'v nr / 


+  °4  i>  + 


5^  f lir  4-c<x  x  ls^o(  %  +^K  c\ 
Z-.  l/W  *(•-**)  1  C'K  Vj-V-J  °J  kC  2(/'^)  °0 


(142) 
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The  stresses  for  the  fiber 


L£iLS(rx-„  \r~  If  ej--ivT-n./v.1i\  q<*  n  )  r  .  r  1 

E-X  11  L —  yhmkW  \  ^  /  -y£)  k  •*._/ 

•*=V*  ■  l<"‘,  3-l 


.r  r /  4-4V>r+71/?tfM  of/t A  ]-r  3_<?k> 

^mir^  ^V'VZ)  \  A-!  +ef(i 


~(w^n 


°W<  l‘Hl-»'JfLJ'‘  *h(!-V 


TjXjJ  J 


_ ( _  ^  r  jo1  .  *  T  /^-r 

^-»I)  V^xr  *e(i->>TJ  ^~° 


t  +  \>z  X 
£I  ^4  " 


V“  v2-  -2+2vx+yL  ,  r  T 

£_  2_  1#/n*L  <f<7->'x;  (  ri  J-U.,+  +(,-*-)  °<k *StJ 

•6,1*  ■■  k-  1.3,  ^ 


^  *&> Ij + 

r  r-jj^fzJl  rj  t  +  °^A  r-  ■}  , 

/^L  ^o-^j  ^-v>x;  X^J 

o(  f  ^  r  4-  ^  7"  i  L  /&<^uo(,  c<m  7 up 

‘»A'L«f().v >r)  ■L'H+I  <f(,-P*)r^-L-YUJ  k 
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t  to-  >i, ^ r-} 

t+v1  x~~  f  r  r-4*^v>-f+n./*i-/)r  / -* ,  r  1  , 

^  *  \£  g,  W  '^?r(  * ;r~'  rJ 


r  r^-^;+Ti/7tf{|7'  -/'<?i;X,v  7'  ]_*. 

L  4(/-v»i7\  Z1  4-(i  - \>s) °^k  At  J 

«&  c“7’^ 

£  ^jiJJ  “"“W + 

I  rt£  v>r  rr 

^i-wrj  ^z,°  '**0-vx)  ° 


w  +•  iiii£  ^ 
°  /-^v>r  r 


£U'_ 


-*g£2ig+ 


nn/c  <K'-sr)  U  *“/  Wnk  <Hi-JT) 


n+< 


uv 

er 


VJ  -  f  f  Is  [J&Lr  +-iirj  1 

**<,*■■■  bv*..L 


p1  r  «jjL.  r  +  t  1  + 

L  4f<  n^;  /-  t>rJ  •"*'• J 


I  *1 


oero,^  ^  oota-Mcf  + 


y 

°*  -n/<  4</-VrJ  A.  J” 

ZT  {fhtL+fc-*1)  *•  +^7r)  A  *  **  40^/° 

k-t,3,  v- 

_  ^r  *  JL/fr  +  _llil£/3r 

i  -  x> x  y7^/o  ^  o  /  ■  <*  v>J  •' 


(143) 
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The  displacements  for  the  resin 


X  /f- 
■^iitiS  *( 


X  /  -  o{  — '  ,TL  T"  \  j. 

4^ ^■v'.+l  4f/-Uir)/L 

£  -,7^1  *0+ 

A*/,*-- 


JT  -f( 

°L 


o/i  f(/-V^5  '  ^  %■ 


f  > 


00  «*  r  jf  -  4  t<^+?i  __  4  7"  i 

=  21  E_  L /7  /  <i-(i-\>*)  ^*1-1  *  'w*/ 


2T  4  -4^-^4-n 


JC  _  4  +  4Va'  4-  71 


x  4-  *, 

tSw  /c^+/ 


^-kJ  c'°°'°k  J  <°^k' 


r 

ZLU 


5^7rJ  -S.+/3^  + 


*  *k\—  T  ^  -rX  °ilr  H 


^S<A 


W  ^/-v^j 


+£&  ^  K  +  <!l  &*,  lo  + 


°<<C  /V1  "1  • 

/^>/Jr  0J/0* 

The  stresses  for  the  resin 


Co  1- 


(144) 


9^'*£  £  k1f(^©-^A/4*><-CJ- 


**V£  /<---/,  J, 


r/  +-4X+II  [ftti )  <4*  A  \r  .  3-29*  ,  -r  1 
+('-**)  W  *(!-»*) JJn+, + 


X 

3w 


f/lii*ffi:Li!L  )//  .5-^y  t- 1 

U  a  /  *(t-va)J,'n-/  y,-**)6* 


x  f(  +-4V*+n(n«\  iA  ],.  _  •3-^vJ  ,  ,j  -| 

jCw.IV  +((.?*)[  n./  +(i-i>B)Jl'n+t  <t-(i  ~Wrylf  + 


A-'M  +  A_  )  // 
<f(  /  -*>u)a*  *0  -**)  J /c?f 


&arzx*7l(j)  ~h 

k:lt3,~ 

•Ct  fejfej*  T‘ '  *>>  k‘+44  4 * 

,  r  w  ^  ^T-2-  -r  _L  /i,t>g  x 

<*»-**  J  +S{,-^)^,0  +i(l-fiT)  ^~a 
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^  •!■]>  K,  ■,  K 


i  jjiji  il 


,  i  i|ji  ii  *jijn  iri  | 


J  Mil U  lU ILL  JMitiPi 


H~VZ 


**V*,'lr=lt3.-L  M 


*  r4-4v*±7t/2i£L\T  _  T  1  + 

(  ^  ^(i-^'a)  J 


2Z-  +.  flirl  w  £/  ] 

-^l/vk'-  <H/-u>:zr)  I  A  / *Hi-\>s)  l(*h-J 


V  r±d^±2Lf7LU)j,  ,  iz*?  v  / >  1  l. 

ef-(<-v,a")  ^  A.  J ^n*i  <?-( #- v,-=r)  '  Wi.  J 


-run- 1) 


,4> 


£{/£  W  Ok  ^  t 


4  cm-°<*}' 


—  x  /  ;  ^  nr  vg_  i^ir,  i+vx  -,x 

Q  A*  2(1-9*)  +  2(t-l>:a)  °  + 


H\>K  xr  ^~- 
E-JT  ^3  “  2_ 


?t -itH- lt^3. 


V~~  f  *  l-<2Vg-ft-  ,  7-  ,  -Z+ZV*-h  ,  T  -f. 

^-fliinlc  4(i-vw) 


jr  u  X  -Z+ZVX-Ti  ,  ty  . 

&**k  <H<-^)^K, ,.,+B^  ^-k^^KhT 


X  -n°(k 


°<nk  fO-^jA.-^4^  tf-O-t 


/&OTL 


o ^  ^  ?l<f) 
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Ac, '  *  = 

jr  r  «//l  T  .1*xJLZ.t  1  . 


x  r_c£i_  .  _  *„»»  ,,  i  + 

•^nr.k  *-  4-ti-)t3)  '•y.-t  /d(i-'/a:) 


JZ  rjQA  , ,  *«”*’■!<'  1  + 

-£->HnU L  +( i -**)  A'**'  *0 •^jr)  nJ 


*  **  J  +&*  /c^lx 
°W  lyJP^k  90-»a)  W 


co<i-  c(^  ^  Cscra,  mcf)  ^ 

44- 

ir  r  y/  _  ^  J^Lr  + 

B^L  *(!-**)  K°J  <HrV*r° 


a1  5^5  /<}  r&r*JL  + 


+  Z-^/3 


(145) 


Equations  (142)  through  (145)  were  then  programmed  in  order  to 
determine  the  remaining  integration  constants  and  then  the  displacement, 
strain,  and  stress  fields. 
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APPENDIX  II 


THREE-DIMENSIONAL  LOAD  TRANSFER  AMONG  FIBERS  IN  A  MATRIX 


Displacements  which  satisfy  symmetry  conditions  (1)  can  be  written 
as: 
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4 

ISJiA  U-4)>+7l)±tih 

+  JoiRttiln  +  'BrtKn)}  MS/l z 

+f{^I„  +  B*rf  L  j  Z^zLWf 


<50 

"I 


u 


^4J  *o 


0-4-p)j-fljtflo  +BaJtKij  /mm-/  ms/ 


S 


+■£0 1  ~  ^  +  +  ^sm 

+  $*»(*-W9)Jt**'  +  flaKA^Z* 


+  B Jf71-' 2* 


flj  (3/]?)  +  A  +  [El  (3-4-p)  +■  BjJ . 


+  fl«  +  Bd  if  r  1 )  f  <54/  < f 


(147) 


60 

+  1 
^70 


4(i~V)l3  —  H  £  ^fl3i^InH  K?H  4  B?/»^  J 

+ B^K*) + (3-#)(y?3^l7f  +I^/f#)}  Mtifti 

I'Wk  I?i  + 

{-/u((#mX,  +BntK,)cfc<(+(/ln& 

v 

+Jik(fLl4 1,-  +  BaM  Ki)  +  /■  (.ft*  •£» +  ^  ^o) 

+ (3-l»)(fl3tI,  +  B  3jKc)}awA  z 

'A{ft*lo+B#(Ko}  *“*/** 

ti b3*J£h]  z 

+  j(i-£v)[fl3+B3MA]z  (148) 
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where  /J^ — /?/,J#  (4^4P~%)  ’'fill  ZL 

B4-n=Bj,m(A-4-m)+B*W  ,  Bt=46-p)B*rB 

and  where  In  and  Kn  are  modified  Bessel  functions  of  the  first  and 
second  kinds  and  order  n  with  argument being  omitted. 

The  first  numerical  subscript  of  any  coefficient  in  the  above  equa¬ 
tions  refers  to  the  term  contributed  by  the  corresponding  Pi  (i=l,2,3) 
and  the  first  two  subscripts  12  or  21  by  Pi  and/or  P2.  The  coefficient 
without  numerical  subscripts  refers  to  the  term  contributed  by  Po. 
Therefore  the  four  Papkovich-Neuber  functions  in  cylindrical  coordinates 
have  to  be  considered  in  the  present  analysis. 


% =f  (M)/  2*  ] 

+  ~/^Kn-i  -[2-8V)JlK» ) 

-  /Ja^[(4'4^xl**»  -  (3~2v)jfa lx] 

-  B5i^[(4-^)x/Cti+>z^+' f 

+ ^z)lJ 
+  5^  [/U  I* + Mj]  ^2  +^(fe«fF+/^J 

-  Z^i^r^y-Wn]}  z  z)  WS7*f 

+  xf|/)^  ^[/Hklo  -  3  A[~pjp  K 

+  ^f-Ic  +  B^K>W^2 

+  {fl^/H^AIr3(l-^)lJ 

-  B^/w[l^K,+30-5)^/{J} 

-  <5)^4  [fljjtlo  +  OOS/iZ  , 

+  [fl*  [#i  -/zi]-BmH¥/z  K‘]  1  Z-^ZJ 
*-  zJ  ] 

+  i[aMj+V,(2-2V)+av^faft  k]Msf + 

+  $f}arVya+M3+B3(z*y+^AJi:) 


(149) 


4o-n 


06  CO 


T  (G  =  £  l0 +  (*-* +*»M X. J 
+ B«^-4+4P)^  /c,r  (*-a  ««)>  K*] 

I*] 


+ ^ »J-^J 


ws/2  ■>-{^[^'I„H+zjrI»] 

?»&-K*h+  ^K*]}  fW 


.  <f  WStf’- +  O-AV^lJ 

+  “4  [Pr1*  ~A  1 1, 


-B, 


'/-■3>?l 


£0S  -ft  2 


+  2x  -'f\frt&-l)jtn~  -  (n+l)J?*~A 

+  f\u»n(-*+M* -B,*nK(*-i)JT* 

f  ftwiifa-i))?  Z2.A +’B3n'?i('M+0)c  z*} 

*  yd^n  7^  ^ 


+  ~i~0~ 


~(BrBjfii#.Ji).4M 


(150) 
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oo 

+  z, 

■krO^ 


=i  i«-  -  m/i j 

+  &*»*  + ( f-«*  J>)/£  Cj 

+  fljl  Ix+r  (I-JP)AX*J 

+  + (b4V)/C  K ’n] 

-/UliU-^n J  +  M£ U+ W 

+  5^-[^3^I,J,  +1$3M  K»1  j  MSpriL 

-lufu’-p y 

rf  4  0-4^1,  +  xO’vXh*  Z* + &  J 

-  %  (i-sp)/I,  +  S/X- 6-^Xj} 

*  <{  Avwtf O0S/C2 

+ f/U^r, -(i-^y^BM[^,^K] 

+£){-/?^-0a’1'*  +  /9j^6*+W*+*)A* 

+  Bin&+l)AX'3  +  Bian(-K-l)(?i-A)Jb~n' 

+  ^*ma-*v+  ^^+/)rM  ^5*f 

+ jj<3  k^k)- 6~3)>)BJm6  f - 

+  5  flal  +  Btf*  +  3Vfl3  -B 3(z*A  -zvffk) (152) 


107 


4-(hVs)^~  «_  eo  f  r  ,  v 

E  V3  ^ ^  $13716 (3-£y~%\tflC-j 

&  sink  I„H  ""  Kv  H 

^(fal*  +fa  KnYfiO-wi^ln+B^M} 
•fafZ  +%r/U2*+&*f 
-£(a(o-^)(^L + fa  K)+fU  i«  +  b^/c] 

•^Wf  -/O-^ifalc+faK j'fas‘f}'4*yM'Z 

+  £^Ti[f\i^)l7" '-  83*  Jf*'']  Z 

(153) 

+  &iz7ia[a^^  K«-rJ^i  l(^]  ^’Bji^jji'^- 
+/%Yl-n+^M^}  ZMv/^cosn'f 

+ 1  [{(fa  (ayil,  -s&z)  -  fa (ayiK,*  f^K) 

*70  *9 

+^(/uj.+bm]^ 

+(fa(Av/IrX^XYBJIY^/Kl + Tj?K*)]f^>i 

+SVJIJ.)  +  Kt  -  aydKo) 
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j 


l 

' 


to  +&/C)  +#(  to+^fC)]^2 
X>  +  1C  j  z^a<C>zJ 

+  f ;  { +<3  A71- 

"l"  ^  [/CC + Bi  +  &t^k\oos^  +'B^Am^} 

+  +  5  to(^3  +83&fJt)  (,54) 

I 

Displacements  and  stresses  as  given  above  can  be  applied  to  any 
particular  domain  of  the  present  problem  by  suitably  specifying  arbi¬ 
trary  constants. 

To  satisfy  the  periodic  condition  of  matrix  and  loaded  fibei;  we 
take  in  (II-l) ,  (II-2),  and  (II-3) 


JV  —  (&  ,  /«3 


IS, 


(155) 


-  T2IfK  -  a- x’%-  TZXf%  n  • 

=  &M  -flj  -0  J  =/,<?  (156) 


Satisfying  the  loading  condition  0 »M't£)=of  ,we  choose  in 
the  expression  for  Q~  (see  equation  II-9)  the  u  so  that  cesuA-n- 
therefore,  '  Z' 


' 
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■tfWiiaUMUiSili 


/I 

>  =ir 


fa  —  l,  3;  ^  jdc  . 


Hence, 


R3M - ^3-4  ' —  0 


c2yX(^-H)  fan  +  —  0 


(157) 


(158) 


In  addition,  the  unloaded  condition  C(*A±L)  °f  domain  II 

leads  to 


/iXZ=J~  4  =  1,3,$-,  jtc.  (159) 

nln = =b£* =#=  0 

tfoi+OflZ, +0-)?)&  =  d 
+avx/£i  =  0 


(160) 


The  displacements  in  domain  III  do  not  have  any  periodic  property* 
in  the  polar  angle  ,  and  we  discard  the  terms  contained  in  the  (J^  . 
and  functions  in  II-l  to  II-9  characterized  by  the  factor  s injttf  * 
or  for  domain  III. 


X  — 
nisnh  ~~ 

tsl 

Brsnt  ~ 

film 

=  BJU  = 

SL 

=-r 

U£ 

P&  ■ 

=  BwS  = 

TO- 

% 

X 

ft  AM  ~ 

11371  = 

=  B*  = 

X 

x 

0 

ft  tot  ~ 

i  — 

(161) 


The  vaniahing  of  the  axial  stress  requires  at  the  end  of  domain 
III  In  equation  (II-9) 


&TV 

^  -At 


&=■  I,  3,  jdc. 


(162) 


TL  vr  TSL  TL 

RsA  =B^  =Ba  =b3  =  0 


flf+Bf  =  o-mf =  o  <l63) 


*As  far  as  a  hexagonal  element  as  shown  in  Fig. 30  (a)  is  concerned, 
displacements  of  the  six  unloaded  fibers  possess  sixfold  symnetry.  How¬ 
ever,  an  elasticity  solution  is  not  valid  for  a  discontinuous  medium. 


Ill 


APPENDIX  III 


Using 
and  general. 


where  V 


and  where 


DISPLACEMENTS  AND  STRESSES  OF  THE  MATRIX 

equilibrium  equations,  Ref.  6 ,  strain-displacement  relations 
ized  Hooke's  law, the  following  is  obtained; 


/  ae 


XV  £ 
xsJ>i  X? 


(164) 


{  ae 


J-V  + 

;ta 


(165) 


_ /  ae 

f~ay*  az 


(166) 


is  the  Laplacian  operator  defined  by 


4-  -LX  4.  X -SlL  4--2X 

x  dX  x*  af  2  a?2- 


is  the  dilatation  satisfying  Laplace's  equation 


(167) 


V2e  s=  o 


(168) 
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Since  (168)  has  been  derived  from  equilibrium  equations  (164) 
through  (166),  we  may  solve  (164),  (166),  and  (168)  for  radial, 
circumferential,  and  axial  displacements,  £  £2>  and  jf 

We  seek  the  solution  of  (168)  in  the  form 


&  — AfoiClnZ 

Substituting  (169)  into  (168)  yields 


(169) 


o  ,  I  dfo 

dp2  p  HF 


(170) 


where  P  -  CA-n^t' 

The  solution  to  (170)  is  given  by 


(171) 


■foW  —  f\\  Ii^p)  +  Bj  K|(p) 

where  Aj^  and  are  arbitrary  constants. 

Introduction  of  (169)  into  (166)  gives 


(172) 


d%  ,  /  a?3  ,  /  g %  ,  l  3%, 
sp2  p  ap  Tp2  3z2 

=  c[fl,  h <rp)  +  B,  K,(f>)]A^L(cJSVnZ  (173> 
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where 


(174) 


c 


t 


The  solution  of  (173)  Is  written  as 


I,  =  f3((>)jm<(cJSOUZ 


(175) 


From  (173)  and  (175),  we  obtain  an  ordinary  differential  equation 
for  with  the  solution 


ts(f)  =  tfff/P)  +/?2l|('p)  +Bz^U6) 

where  A 2  and  B2  are  arbitrary  constants  and  where 


The  function  in  (177)  is  given  by 

M)  =  K,(fil (H8) 

in  which  primes  denote  differentiation  with  respect  to  (0 
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By  the  definition  of  the  dilatation  g  ,  we  have 


Substitution  of  (169),  (175)  and  (179)  into  (164)  leads  to 


+  Wp)  +/U  ^  2 

(180) 

1 

»h.re  -f/f.)  =  Cl,(f)  +^pfl|<'f)+0/«f,(p)] 

(181) 

Vp)  ■ =  c  K/(p)  +fp  &G  (?) + 

Assume  the  solution  of  (180)  in  the  form 

£|  =  ^(fiAitP^AwOlnZ.  (182) 

Introducing  (182)  into  (180)  and  changing  the  variable  result  in 


(e4;)"+fkp  <hf-0+-p)(H) 

—  fli-f/f)  +  8cf+((>)f /?z^^  <183> 


115 


i 
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The  function  Q0  in  (186)  is  an  arbitrary  function  to  be  deter¬ 
mined.  Now  that  the  $2—  plane  is  the  symmetrical  plane,  we  have 


which  requires 

(j.CfMO  =  0 


(188) 


Stresses  ,(£•■>  >  °f  the  matrix  are  related  to  strain  £ .  >  by 
generalized  Hooke's  law,  as  follows:  4 


*(I+)T) 


■gj/e  7 


(189) 


where  1®  the  Kronecker  delta. 

By  using  displacements  (175),  (182),  and  (186)  together  with  strain- 
displacement  relations,  stresses  (189)  can  be  written  in  the  form 


07;  —  T;;  (f)Ato'<(AwO(*  2. 

(Tiz.  —  AAn&inZ 

(J73  = 

<J^  =  %;([)<**<{  CJSOtnZ 
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where  /T^j  ere  functions  of  p  only. 

It  is  seen  from  (175),  (182),  (186)  and  (190)  thnt  displacement  com¬ 
ponents  £  and  and  stresses  (J^  ,  ,  (J^j  ,  and  {J73  characterized 

by  a  factoS:.***^  vanish  at  the  %2  —  plane  (or  the  neutral  plane),  are  sym¬ 
metrical  with  respect  to  the  plane  (or  the  symmetrical  plane),  and 

have  maximum  values  at  the  >Jf2-pl*ne  for  constant  X  and  2  .  The  dis¬ 
placement  component  £2  and  shearing  stresses  0/2  and  characterized 
by  a  factor  cotf  are  symmetrical  with  respect  to  the  ^2 -plane,  vanish  at 
the  'Jfjr- plane,  and  reach  their  maxima  at  the  #2-plan'i  for  constant  X 
and  2  •  Obviously,  these  properties  satisfy  the  physical  requirements  of 

the  bending  of  the  matrix  in  the  present  analysis. 


The  arbitrary  constants  and  B.  (1*1,2, 3)  in  the  expressions  for 
displacements  and  stresses  are  to  be  determined  by  boundary  conditions. 
The  outer  surface  of  the  matrix  is  assumed  to  be  free  from  stresses, 
and  hence  we  have 


07;  (fi/^f/ —  0  cm) 


where 


(192) 


The  continuity  condition  of  the  interface  on  the  components  of 
displacement  in  the  directions  of  Cartesian  coordinates  ^  ,  and 

2  requires  c 


z)am<{  +  z)wsf  =  U.  aUoi^z  (193) 

2)  =  -  av,  atn  MS  O',  H  /Un/(^ 


where  P5  —  0kg.  (L 


(194) 
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-liitfjfain  hi  * 


a— ^ 


The  condition  (193c)  has  been  obtained  by  assuming  that  the  angle  of 
rotation  of  a  fiber  element  is  small  compared  to  the  unity  and  that  a 
normal  cross  section  of  the  fiber  remains  normal  during  bending. 

Substituting  (175),  (182),  (186),  and  (190  a,  d,  e)  into  conditions 
(191)  and  (193)  yields 


/T/;^P/)=0  i  =  b*,3> 

fi  (ft,)  =  Vc  2=1,3.  cm) 

'&(&)  =  -**£ 

The  constants  A  and  B  (i-l,2,3)  in  displacements  (175),  (182), 
(186),  and  stresses  tl90)  cAn  be  determined  by  solving  equations  (195) 
simultaneously.  Each  of  these  constants  is  then  characterized  by  the 
presence  of  the  amplitude  -tT#  as  a  factor.  Therefore,  stresses  at  the 
Interface  can  be  written  as 


Cm  (ft  Z)  —  /Vff  (Til  A+rt/^ 

Z)  =  U  (J},'  OflStf  A+rt,0(n2  (19( 

(J73  Z)-  U  &  A^  OOS  0/n  Z 

where  0//  ,  OT2  «  end  07,  4re  function,  of  ")?*  ,  0^,  c,  b,  and 

Ci  yt 
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APPENDIX  IV 
ENERGY  SOLUTIONS 


A. 


ALTERNATIVE  WAY  TO  FIND  SMALLEST  BUCKLING  LOAT)  AMn 
WAVELENGTH  OF  FIBER  FOR  MULTIFIBER  CASK  “ — 


BUCKLING 


Hie  strain  energy  due  to  longitudinal  shear  along  the  fiber 
surface  can  be  also  found  as  follows: 


~  I  07*  e*  dv  = 

zjv  11  i- 3 


3rd 

■ZGr 


fL  I  f  l  ,  _  Si thtSc<,4)  \ 

t  K  a  \  +s,4m)  5m  * (  ^ 

R’/ra)  ( /  -  cosA^-Wi), 

"  (ioU)  +  s,„J, u*  J)  K4)5i »J,K A) 


_  ^S.V/ra)  si»l>u«.J) .  .. 
■  CBS  (<>(„*)  A g,  Ja. 


12- 


:  AtlS,.)  4  (M 


2-<*rA  +  S/fc  li  ^  /f 

_  i?  *LJ 

+ *  I 

~  f2^/+  (E^i3  _  -f  i)  $  ihli(244) 

-  jr(<M)cah(M,ty  c&bh 
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+  ± 


Si  nil 

^  24rt4  ■f 


2.  f 

)  *  f*)si'nliU<*tty) 


-  g  ((**«*$  1 3-)  S'ml)(j.o(tta) 


■f  *f 


J2*j  ( I  -  C*sh  (24 J))  ft 

(MA  +  smk^M'l2 

-  (^»4)  —j^Ccsl)  (2o(n\)) 

ijf  emit (^„ a)  -Z(o <hlpffJt(e(„«t)l'\ 

_  Zd„<{  Si '^(zttU)  r,  ,  .,  . 

fcwj  +s;^j)f[^  <«t>c*kf*w 

-  ^(24, a)  -  -Jj:  Si»!)f2c*„L) 


-L  1.  ( I  C't%h(-2o(hm);oiyin\4ai.d)  \,.  .a  *  .  .  .. 

f  s  J  WW 

-((^a)l-fz)c^(2^) 

- (24 -j-Jffon0)  Si*h Ue/h ^Ojj- 

(197) 
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. . minus 


P' 


Then  the  expression  for  buckling  load  of  the  fiber  is 


—2 


?;  -  e’i* + ■  r 


ZGTti  a 


(198) 


where 


} 


of  equation  (197) 


Combinations  of  equations  (197)  and  (198)  yields 

fcr  —  octE jx  +  I  s/4 

6 U3(2Kj )-8<xJ+ SotJceAfrJl) 

—  4  fad)ce3k(j2c(Ji)i2Si'4  (2**4) 

+  £Vf  f ?  j 

t  £-2  S  J)  -  8<*d  +84 J  CtohfaJ) 

-4  slJ,  U*d)  <^texJ)+2$UU*J) 

-f  2.  Sinh{Jdd)Ct$l)  ^j. 


(199) 


SMALLEST  BUCKLING  LOAD  AND  BUCKLING  WAVELENGTH  OF  THE  FIBER 

The  smallest  buckling  load  of  the  fiber  can  be  obtained  by  minimizing 

Pi  of  equation  (198)  with  respect  to  a  •  3P1  /da  =0  goes, 

cr  n  cr  n 
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zct'„  E*!*  +  E*irr«.  [ [(2 Jf-  ^  X«t,(z*J)cU(z*J) 

4  £>4»  Z  Sink*  (2o(h4)~^ 

*  £ -£SihL>3(2 o( „  J)  -8o/hA-f~3*«<l C&h faj) 
-4  51  h/i (-2  <*„<{)  <ushfad y+zsinlt  fad) 

4  2  S I  hi  (2o(b4)  CUsi)  (2»in  *r] 

-  £ \<l  j [  {2&T*k  -<**' '  sUl^U^J) J- 

,j^-~2.5ihJ,JU*„J)-&aS«‘i!  +8ot„Jl  C*sh(£.o<nJ) 

-  4  51  */,  U^4)  &*.<()  -h2SlH.il  U*H  j) 

4  2  Sink  (2o/d)  &rshlUo/„J)  J  0 

•  ^-/lJs\dLUo<oJ)  CGi,(j>o<J)-8((-h2<l  tc*>kUo(J) 

H^J^hSikli  -Sd!  sinltLfad)-&4  t&hfuj) 
4f d  c&hfaj)  +4  4  e&kfaj) i $4 sihlffad) avfa'Uj ! 

f 

^  EXJT^  • 

*j^~2Sl  hJj  (j2.o(d)  ~~8^hd  ~h8o/^4 Ctf&k  (2e<d) 
-^sikh  (2^i)cjs^faJ)  i^S/ni^d) 
-h2S\h.h(2a'd)  ^^(2o(hJ)J 
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"  4  f+2)sihl)  ~^ctJ>Cch  h  C&J) 

f  j-WC&U^a)  -  -f 


+2  C^li  • 

*jj-(  H  i1)  5i4^hU  +  ]jf/ 

(WZ)  3m/»  (^y^M-i-z)  odi(M,$ 
-j-  ccsk(j&(hl>) 

C&l)Uo(n4)  -f-o/haz5ihlf(^^) 
z 

-°/Jh3  +  o^a3  j 
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±  _  (**£  1 
*  3  3 

)  31k  l)  (2dhjy) 

f  f*)  £&hfahl))- 

~  jr  (fd^Sihl,  (id*  ^-jr{(2d*(^ft^(u^ 
-j-  OssUld^)  -d^sij \(id*[>) 

&cn>l)(2d„<x)  -j-dh0LSinl)(2d*a.) 

+^HZt3  -<***?  J 

'  UotJ, )  -fax)  Slul)  (Z4hh) 


—  \Zot*  U>~  C«i  U*d))  * 

*  o(nl)  Oftl,  (2<dH i)  -  dh  aZ  C&h  [Zo(*ck) 

-~/j- (dub  -f  J 

~\zJi  J)  -j.  Jfa et*ash&J)l 
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•  <^3li  ( 

—  jj;  Sihi)  (2o(n  h)  *f  jj-  5lnli(-2rf„a)J 

^ )J ' i ^ / H t (2*h b) 

-ot^S  Uh(2<4)J 

-j~g\zc)  CmhC2o(hJ)(l-CC5li  (24*4)) 


2&  Slut]  (jLoCh  4)7# 

'ff('2^nl))X'f~  2)a%lt  t>) 

-((. 2ot»a)l-j-  2)Cj2> \(2  of^j-4  <  h  Si  hi  fait) 
-f-  Jf  o^M  £LS«  hk  (SZo(n  a)J 

f  g~  (l~~  C^&h  (^°^n  d))  SinJi  {-2dhd)J» 

•  Cal, 

* 

~^d*(K  Cc&U^) 4(^ V t)sMfa$ 
-41  to„t)  -gosJ^Cs&hu^L) 
44wSi4  £-&*„«)  +  $<  <?  cali(Mh$  J 

+  EfiW  u+p5  I  z[(uJ)^(Mjj-\z\  -z\  tell&Jjj'o/''. 
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•  --2.5/tti  -f  %o(hJ  C^h(li/hJ) 

\ 

-4sinh  (2«J)  &sh(*<Jj  (2**4) 

-j-jz&i'nl]  (£*hd)cc3ljZ(MHdf )J 

—  d)  —  Si  nh(2^  J )j  .  o(H  ^  , 

”-2  Sink  (<2o(hd)  ~8°^h^  ~t  Os>l)(2o<nJ) 

s 

Sin!)  (c2<xd)c45lf(t&d)  j-JZSltth  (gXhd) 

“M  Si  ah  (jbsJ)  CGhlCfokJ)J 

~~ 2|  d)  —  sikIj  (c^^nd)J  •  oc\  i  • 

*| ’’’"•^  SihJi  (2°(^  d) ~~ S^c/^d  j-  8°/h<\  £&h(jfo(h<l J 
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L 


4  g-  Q—  Cg*>Ii  d  ))  3ihIi  (£o(»d)  . 


(fa"  bJi^Ozli 

i  t>) 


o 


(200) 


From  equation  (200),  we  find  the  buckling  wavelength  ^  =  2rr  /«.  , 

J 

Therefore,  the  smallest  buckling  load  of  the  fiber,  (P1  )  .  ,  can  be 

obtained  by  substituting  lcr  into  equation  (199).  min 
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B.  BUCKLING  OF  A  SINGLE  FIBER  IN  A  FINITE  MATRI)  UNDER 
AXIAL  COMPRESSION 

INTEGRATION  CONSTANTS  AND  LATERAL  DEFLECTION 

By  applying  boundary  conditions  (111)  to  equations  (83)  and  (84), 
and  solving  them,  we  obtain  the  results  of  integration  constants 
as  follows: 


A*- 


-  4*2- 

7TOJ 


“(fe 


C, 


D  = 


1  rW)  ^  -  /V  | 

d yi  ~~  Si  hJj  d n  J 

!  -f  Siyili  dn  C45I1  d h,  j 

i  $1  ftii  dn. 

tj 

1 

v7lV  L  dt  -  SiXJ, 


where 


=  Ipa-a.) 


(201) 

(202) 

(203) 

(204) 


(205) 


The  relationships  between  displacements  and  stresses  of  two-dimensional 
elasticity  problems  are  quoted  as 


of 

vt  - 

E" 

*§T_ 

sf  - 

dz 

Ex 

E* 

(206) 


(207) 
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STRAIN  ENERGY  DUE  TO  LONGITUDINAL  SHEAR  AT  THE  INTERFACE 

The  strain  energy  due  to  longitudinal  shear  at  the  fiber  surface 
is  defined  as 


•  os' Az  A.  (2ii) 


Integration  of  equation  (211)  gives 
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BUCKLING  LOAD  OF  THE  FIBER 

Total  potential  energy  of  the  fiber  during  loading  is  defined  as 


Te  =  Ut  +0e  t  Op  -t-Oj-W 


(213) 


Combination  of  equations  (70),  (208),  (209),  and  (212)  yields 


'7t>i 


where 


X 


4 

r  i  / 


J  !  l_\  r.  b-x'i- 


■XxK  a1  V  K 


UTJhX 


'Tpj 


hi 


of  equation  (212) 


Minimization  of  total  potential  energy  with  respect  to  5^n 
gives  critical  buckling  load  of  the  fiber. 


(214) 


(215) 
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I  or 


*  +  x 


^(tTc  a.  v°^v 


(216) 


where 


H 


—  'TUL 


(217) 


L 


Substitution  of  equations  (208),  (209),  and  (212)  into  (216)  results. 


Pc 


X 

cr 


#|^2.(e/i1  —S)kiid^)(02>kdk'~  Oj 
•  (  dp,  -  S)hI)(Jh)  C  (^3^  dp,~  oj 


(218) 


SMALLEST  BUCKLING  LOAD  AND  BUCKLING  WAVELENGTH  OF  THE  FIBER 


The  buckling  wavelength  of  the  fiber  is  determined  by  minimizing  P 


cr 


BP 


with  respect  to  ;  i.e.,  °  cr  _  q 

Ban 


From  equations  (212),  (215),  and  (219),  we  have 


(219) 
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c.  BUCKLING  OF  A  SINGLE  FIBER  IN  AN  INFINITE  MATRIX  UNDER 
AXIAL  COMPRESSION 

INTERFACIAL  FORCE  PER  UNIT  LENGTH  PER  UNIT  LATERAL  DEFLECTION,  k 


Now  let  us  take  b - or  the  matrix  is  infinite  in  size. 

Thus,  it  is  more  convenient  that  the  solution  of  differential 
equation  (79)  takes  the  following  form: 


$  =Z^z[Ain  e^-f 

+  (221) 

The  stresses  are 

£>o  f- 

C  =  l[(  6,  J  1 6,^)^ 

'Sin  (ol^Zj)  (222) 

qT=~2[A„e"  fA.Me**' 

'7>€L\ 

+  CM  €  ^  +J)Jo(«tye  ^Jsinfaz) 

(223) 
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-  J  ((A,.+  eje"+B„Me^‘- 

(22M 


Since  when  r  -*  =°  ,  all  stresses  have  to  be  finite,  then 


Am 

&IK  =  0 

and 

C,H  = 

Dm=  0 

Therefore,  the  displacement  in  lateral  direction  is  found  as  follows; 


(225) 

(226) 

(227) 

(228) 
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(230) 


or 


nab* 

t-P* 


I  ^ 

->3^1 


(231) 


SMALLEST  BUCKLING  LOAD  AND  BUCKLING  WAVELENGTH  OF  THE  FIBER 


The  expression  for  buckling  load  of  the  fiber  is  derived  as  follows: 


Pi 


°CEPr‘  4«.-' 


(/  -^r 


(232) 


In  deriving  this  equation,  the  value  of  k  in  equation  (231)  has  been 
introduced. 

I 

Again,  a  is  determined  by  minimization  of  Pcr  with  respect  to 
a n  of  (232);  i.e., 


(233) 


(234) 


Substitution  of  (234)  into  (232)  gives  the  smallest  buckling  load 
of  the  fiber. 


which  yields 

N  -2E 
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BUCKLING  OF  A  MULTIFIBER  COMPOSITE  DUE  TO  RESIN  SHRINKAGE 

Because  of  the  smallness  of  the  cell  diameter  compared  to  the 
length  of  the  fiber,  we  can  make  some  simplifying  assumptions 
such  as: 

1.  Fiber  cross  section  deformations  due  to  stresses  perpen¬ 
dicular  to  the  axis  of  symmetry  are  neglected. 

2.  Cross  sections  remain  plane  after  deformation. 

3.  The  shape  of  a  cross  section  is  maintained  after 
deforma  tion . 

4.  Cross  sections  orthogonal  to  the  axis  of  symmetry  remain 
orthogonal  to  the  deformed  axis  of  symmetry  (no  defor¬ 
mation  due  to  shear). 

Figure  64  shows  the  deflected  axis  of  symmetry  of  a  fiber  and  its 
associated  hexagonal  cell  in  the  z-direction.  An  arbitrary  point 
P  at  the  nondeflected  axis  has  the  distance  z  from  the  origin. 
After  deformation,  Point  P  rests  at  P  with  the  coordinates  £3 
and  Similar  is  Q«  the  end  position  of  Q  after  deflection.  The 

point  X  is  an  off-axis  point  of  the  elastic  medium.  It  will  be 
at  X  after  deflection.  Because  of  assumption  4,  the  line  XP  must 
be  orthogonal  to  the  tangenial  line  of  the  deflected  line  of 
symmetry  at  P. 


J-fciLj, 

V  « 


<1,  My 


Figure  64.  Deflection  Geometry  of  a  Fiber 
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The  strain  of  the  element  is 

As, 


fM  - 


xy  ~xy 
"xy 


_  ;  r 


A* 


or  for  the  limit  Az  — $•  0  , 


ds,  /,  X 


=  7¥V'  J 


- -/ 


(238) 


For  x  =  0,we  get  the  strain  of  the  line  element  along  the  axis  of 
symmetry 


and  therefore  (238), 


£3(x)=  <?,(<>) -x 


(239) 


d_s  I 

~ds  ~f 


p  is  the  radius  of  curvature  of  the  deflected  axis  at 
Because  dtp  =  ,  it  follows  further  that 


P. 


SM  =  ?,(o)  -  x 


d 


t 


d& 


(240) 


e^x)  is  only  a  function  of  x.  Therefore,  all  fiber  elements  with 
the  distance  x  from  the  plane  normal  to  the  plane  of  symmetry 
(in  Figure  6'+  the  plane  through  line  PQ  normal  to  the  paper  plane) 
experience  the  same  strain.  The  strain  in  the  x-  and  y- 
direction  is  caused  by  transversal  contraction.  Quantity  v  is 
the  Poisson's  ratio  of  transversal  contraction  of  the  material. 


We  get 

£Wv)  *  ~V  <53W 


(241) 

(242) 


Considering  infinitesimal  distortion,  we  can  write  the  strain  as 
follows  (putting  63(0)  =  6) : 


By  neglecting  the  shear,  we  obtain 


V  -  D 

23  '3 


(243) 

(244) 

(245) 

(246) 


The  variation  of  energy  is 


Putting 


(247) 


(248) 
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Performing  the  integration  over  x  and  y  ,  we  get 


where  A^  and  A-  are  the  cross-sectional  area  of  the  fiber  and 
resin  in  the  cell. 

Fur thermore, 


r 


( 


0 


the  areal  moments  of  inertia  with  respect  to  the  neutral  (non- 
s trained)  plane. 

Because  of  symmetry 


f  X  <4*  <hj  0 


then  (250)  and  (251)  become: 


0 


The  equilibrium  equation  in  energy  form  for  a  fiber  resin  cell  will 
be  written,  after  introducing  (258)  and  (259),  as  follows: 


(258) 


(259) 


(260) 


Due  to  deflection  of  the  axis  of  symmetry,  the  projection  of  PQ  on 
the  undeflected  axis  of  symmetry,  the  length  A§,  will  be  different 
from  Az,  the  original  projection  of  PQ.  The  difference  A§  -  Az  means 
the  longitudinal  change  of  length,  Aw,  of  the  axis  element  Az. 


AU/  -  —A*  = 

3 


Taking  the  limit  Az  -»  0.  then 


(261) 


Because  of  (239), 

^  is. 

ds  ds+L 


cL  s 


C*n  {j>  "=■  (_  I  ~t  E)  '■f 


and  combining  the  foregoing  equation  with  equation  (261)  gives 


thvJ  ~  j (/+■£)  cat>Y~ 


(262) 


The  total  change  of  length  will  therefore  be 

a 


w) 


,  I  j(i  +  tj  f  •"  /j*  ^ 


(263) 


'0 


We  assume  that  the  fiber  is  restrained  to  the  matrix  such  that  the  longi¬ 
tudinal  change  of  the  fiber  is  equal  to  the  longitudinal  change  of 
the  matrix.  Therefore  it  will  be_J~ 


and 


(264) 


L 


( /+  a  a*  f 


-(/+/f 


(265) 


The  solutions 


£  -  £<W 
Y  ~  f&j 


(266) 

(267) 


of  (260)  have  to  satisfy  equation  (265), 
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We  consider  the  case  =  0  and  get  the  following  equations: 


with  the  restraining  conditions  (265), 


The  Euler-Lagrange  equations  of  this  variational  problem  with  two 
independent  variables  are 


+  7  -=>0  (269) 

6V+  (270) 

where  we  use  the  abbreviations 


<p=  EIAT+BTAT 
<2  =  EtIt  AEJrIJ 

m  -  -  £rAI[iF  >  o 


and  the  parameter  X  being  the  Lagrangian  multiplier.  Solutions  of 
the  set  of  equations  (269)  and  (270)  must  satisfy  (265). 
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From  (265),  we  can  determine  X  by  substitution  of  e  from  (269)  in  (265) 


-A  ^  a?  =(  l+(i)A 


(T,-!£)  J I?* 

- _ (274) 

/ 

0 

In  order  to  give  a  physical  interpretation  of  X ,  we  consider  the 
resin  around  the  fiber  not  shrunk  (3-^  =  0).  Instead, we  apply  an 
axial  force  =  F  acting  at  z  =  l.  Then  we  have  to  apply  (260) 

with  0H  -  0  Land 


-£  F^it-  UF  ) 

AA-^t  A  *7) 


The  variational  problem  for  this  case  will  be 


-  0 


from  which  we  get  the  Euler-Lagrange  equation  of  the  form  (269), 

(270)  but  X  =  +F.  The  parameter  X  therefore  represents  the  axial 
load  of  the  fiber  which  is  generated  by  the  resin  due  to  shrinkage. 

We  can  assume  X  to  be  known  and  substitute  e  in  (270)  by  e  from  (269). 
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The  following  differential  equation  in  P  is  obtained: 


After  multiplication  by  dP,  we  can  integrate  (275)  at  once,  ard 
obtain 

i  -  :L  \  f  ,  '12.  }  .  ,  x-  _  O  ,,  _  q 


4  fp  Cin2f 


(275) 


(276) 


We  assume  P*-  =  0  at  p  =  Pmax1;  therefore, 


C  -  -\(i 


2-  )  Or>  </  f  2  f 

?)  W  4P 


(277) 


and  (276): 
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The  integral  in  (279)  cannot  be  expressed  by  elementary  functions. 
Using  the  abbreviation 


the  integral  to  be  determined  is  as  follows: 


(284) 


(285) 


(286) 
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j  jd-^{ 


Finally,  we  substitute 


The  integral  in  (291)  is  an  elliptic  integral  of  the  first  kind. 
Its  inverse  functions  are  Jacobi's  elliptic  functions.  From  (291) 
we  get 


d  $ 

- 


(293) 


/7VM  ^  /Jvv 


(*  -  C- 


a 


-lit? 


Re-substitution  gives 


UJ  — 


i 

(/•t-2^ 

/ff f  i 

r 

2  « J  ^ 


f  -f 

/»  2y 


(294) 


Solving  for 


/?Vvi  gives 


/7^2JL 
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/rvvi 


( /+ 


h  /° 
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z'yw'  j£ 
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±  /jv*p-<:,)x 


Z±H 


‘y*  -  ?*~f 


.  (295) 


with 


K-\ 

(f  -1/ 

Because  of 

/  ‘f  ^  »  ' 

-  Z  Cy}, 

Vz  -  £ 

follows  (295) 

^  i  />w\  y 

(296) 


-  jpf /f  _ 

(^K  y 


(297) 
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We  obtain  function  §^(z)  by  quadrature: 


q  y  , 

I  —  — - —  /Tw  y)  ttyr 


sinPP  and  costy5  are  known  functions  of  z  according  to  (301). 

We  are  now  also  able  to  calculate  X  according  to  (274).  The 
calculation  of  §1  and  X  depends  on  the  calculation  of  the  following 
four  definite  integrals: 


=  f  C^>Y  ^  ~  ~  2- 

Jd 


J= 


-  2  f/*~\  ~  on  J 


*  A  C 

-  f on  Hy  *  8-8  f 8  Mw 1^-  ^ 


We  have  to  solve  the  following  four  basic  integrals: 


(310) 


(311) 


/ 


X 

J_=  f 

IV  J.  i 


cp 

For  that  purpose,  we  have  to  evaluate  cos-^  from  (295).  It  is 


(312) 


(313) 


hZ<Y  -  + 

^  l  —  y  />w2|zx  +  kj 


I  +  2y  —  y /^z jjrK  +  kcj 
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Because  of  (292)  and 


follows 


(314) 


All  integrals  J^.,  J.^,  can  f°un^  in  a  closed  form. 


It  is 


Q 
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e 


For  the  first  integral,  we  substitute 


and  get 


f  <*C 

J  1  -  r% 

K  7 
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1}  (  ;  k,  jfjcy  j 


no 


L 

K^jk/n)  J '  (i  i-»ix)\[(i-t'-)n-kti*) 


is  the  elliptic  integral  of  the  third  kind. 
The  integral  (319)  therefore  is 


£rTH*f -rffc 


Next, we  treat  the  integral  J 


II* 


Here  we  substitute  cn  (  z  +  K  )  =  t 

=-*  />w(xitK)  ^  ^ 

Av^XS  +  K)  =■  I-  <X?(V^+^  = 


We  get 


cn(n£+K) 


'y 


\ffo^)  tucJxui.  (cn  (K+k£)\I 


The  third  integral  will  be  solved  by  using  the  transformation 


(‘H  z  +■  k")  —  t 

dt  =•  96  \lO-i2)(l-  wj  dz 

Sn(HP+  !\) 


(323) 


Since  the  rational  factor  of  the  elliptic  radical  is  of  even  power 
in  t,  we  can  rewrite  (323)  to  bring  into  a  more  common  form. 

But  first  we  have  to  expand  (323)  into  partial  fractions.  Setting 
for  abbreviation 


Became  of 
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follows  (325) 


s*y(K-rx-£) 


S*  (Ki  nJ!) 


-  z 


SnK 


SH  (K  +  K^) 


f(i-»%iV\/<)-t*Xi-M)  }  =  (P^  \l 

&h*  K  J 

f n(s*(K<-**)j  k)  ~  k)] 


_  a.  fn (in (*+■«■*)/  -*> k)  ■ 

^  l  shCh'+xS) 

±  f  At _ 


(326) 


In  order  to  evaluate  the  last  integral,  we  shall  derive  a  recursion 
formula  for  this  kind  of  elliptic  integral,  since  all  elliptic 
integrals  can  be  reduced  to  the  three  basic  kinds  of  elliptic 
integrals  and  rational,  logarithmic,  and  inverse  trigonometric 
functions  in  t  and  the  elliptic  radical 


16  8 


First  we  expand  the  polynomial 


(l-^)O-k'L')  =  l~(H-kyX--bkt-i\ 327) 

(l-n'-i1) 


in  the  power  of 
Setting  t2  =  §,  then 
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*  A-J>  tAiJ-'  *  A‘x" 


(335) 


The  left-hand  side  integral  we  integrate  by  part 


a 


oib 


ey 


_ u 

y  J1)  (i-n'ty*1  y 


f 


Evaluating 


l  M 

(336) 


dX~ 


ctt 


»  (AY‘ 


and  then  (336)  becomes 


=  -* 


»>*<*■ 


Na- 


f  U  cUt 

Jo- 

+A,o-^v)  f  a,(/-^T] 

-  A,  (!-»'?)  -  2At(l-»%*t/J  at 
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or 


r  «  cLt 


d.t 


ZvA, 


+  A,(2v-i)(i-^)  4-  2(/-l)  A2  (l-n't1)  ■ 


■h  f _ 

JO-”H7y 


f2vA*  +  A,  (2v-  /)(/- 


—  a 2v A«  3'v+i'4‘  (2\fAff  A/ (2v-ijj  J, 

+  (A,(2v-l)  -  2(v-i)A2)  Jv_,  +(2v-2)Az  Jv.2 

(337) 

Replacing  the  left-hand  side  of  (335)  by  (337),  we  get  the  following 
recursion  formula: 

2v  Ac  J "vti  +  [(l  2v)A0  +  A,  (2v—/)J  Jv 

L  2(v-I)A,  +2(v-i)A  2.]  Jv-i  ~  (2v- 3)  J"v-z 

=  _ 

(’/-  in'i'-y  (338) 


Because  of  (329),  (330)  and  (331)  are 


2  V 


2m Ac  =  (n4 


(339) 


(2v-l)(A,~  Ae)  "  ~^T  (”*  +  3k1-) 


2(v-f)(  A;-A,)= 


2(v-l) 


(340) 


(3>k‘l  ~  ^ 0  +  ^))  (3*41) 


“(*'-*>  A,  =  -  ^ 


2v-  3 


/f 


(342) 


From  (338^  we  get  our  integral  in  (326)  immediately  by  setting  v  =  1; 


2A,Tz  -  +  (A0  -  A i)  Jt  A 2  vJl / 


or 


J,  - 


1 


t 


2(h‘  htt1')  Ao 

We  have  to  modify 


-h 


Ap-A/ 


T  -  Al-  J 

2  A-  -/ 


2  f\o  2A0  ■  z‘^°(343) 


JL,  -  /  Xirrjl  <st.i 


into  the  canonical  form.  We  get 


M  -  A 


1/ 


ot£ 


3 


h 


K 1 


•? 


(344) 
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E  (t,  K)  is  the  elliptic  integral  of  the  second  kind  Now  we  know  that  the 
elliptic  integral  of  the  third  kind  is 


=  I~T (t;  ~n "  k) 


therefore , 


(345) 


From  (329),  (330),  and  (331) 

It  l<T) 

K  v _ 


and  (345),  finally 


nAt  va-t'W-k'i ') 


0-rritT'J(l-t'-)0-kvl'y  2(I-hvi +kl-»'-(l+k%) 


2(>n«bk''-n Y/^A'1'))  '  ' 


- *1 _  £  (iik) 

(h*  f  k'm-»'’0+k%))  '  ' 


K^bk1-  *x(i+kx)) 


(346) 


Finally, we  are  able  to  write  the  solution  for  JTTT  from  (326)  with 
(396)  as  follows:  111 


-r  /  0  /  +  <j£.\2  f  /)  c,n(Ktnl )  cLm  (K  -t-xt ) 

JE  ~  [P  f  I  j  *  + 


- 2r\x(\+kx)  r  n/  v 

7*7*  1 


2'K(h‘1bk'-’-n'~(n-k>')) 


[£ (sn(K-ml)j  k) 


o 


After  this  rather  extensive  evaluation  of  the  four  integrals,  we  are 
now  able  to  give  the  analytical  expression  of  the  deflection  coordi¬ 
nate  T)  of  the  deflected  line  of  symmetry.  According  to  (305),  (307), 
(309),  (311),  (313), 


^  =2( ~  jr  [2  J^(z)  -  4  Js (z)J 

=  2(1  -b  ^  (?)  (350) 


We  can  now  express  and  JIV  explicitly  by  (322)  and  (349)  and  get 


('/t? 

P  ^  J  V 


c*t.(K  +  *% 


?  f 
_  2-2  £. 


>) 


PK  V  f  ’s'  l-H%  -h  C,nl(K  +  ?i r2  ) 


or 


cm>U~  (Vr+e*  (*+**)) 


\ 
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- 1 


4"  «r 


cn  (K  +  m  ) 


I  ~  h*  -h  env(K+'>x:z) 

(351) 


(351)  represents  the  deflection  curve.  The  periodicity  of  this  curve 
is  the  same  as  the  periodicity  >  (as  found  for  cp  (z)  ),in  equation  (302). 

We  shall  introduce  the  original  physical  parameters 


-  <b 


CV~D  i 


l-h"  -  ^  or: 

*  **" 

EL  -  il  -  I  -  ,  - 

y  -  2jt  1  ~  a  ~  '  - 


P-R-2 


I  - 


je+p 


o 


(352) 


(353) 


(354) 


I  -h  &  =  - - - - —  -  -■  ■  ,  —  —  f - i 

2(P-R-a) 


I  +2a  = 


P-P-  A 


(355) 


(356) 


\j2PQ  \[r~* 


We  have  to  determine  cp  and  \ .  Cp  can  be  determined  from  (304) 
and  (262):  max  max 


(357) 


and  because  of  the  assumption 


•A 


Then,  \  is  determined  by  means  of  (274),  (306),  (308),  (308),  (321), 
and  (347) . 


(P-R)[/- 

5 '  i-  -  f 


(360) 


In  (360), the  expressions  (321)  and  (347)  apply.  (360)  is  a  transcen¬ 
dental  equation  in  X,  since  J^.  (1)  and  Jj__  W  are  functions  of  X. 
In  order  to  evaluate  X  numerically,  an  iterative  method  should  be 
applied.  In  order  to  do  this,  we  have  to  approximate  the  functions 
J_  (£)  and  Jj-j-j  (t) ,  so  that  the  first  value  of  X  becomes  reason¬ 
ably  accurate,  but  this  means  that  (321)  and  (347)  must  be  specially 
investigated . 
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The  complete  Integral  of  the  third  kind  tt  (  -  n  ,  k  )  can  be 
expressed  by  K  (x^k)  and  E  (x^k)  as  follows: 


rr(-»l,k)  =  K(k) +  ^==^=  [E  ( fc;  k)  K(k) 


-F(£;k)EM) ] 


(365) 


for  /  n  /  <  /  k  /,  which  is  the  case  in  this  problem.  For  1> 

Jz  -  fp  +  ?)  (^==~  ==  [ E(£;k)K(k ) 

-  F(^;k)E  (*)] 


(366) 


The  integral  (347)  is 


X  -  Jlf-lv. *■■))«{  (4N-‘>  KM-*'*”'-  7kV 

-  O'**7"'-  [WMW 


-F(%;k)  E(k)]  -(tN-l)»'-E(k)J 


(367) 


Therefore,  (360) 


A  - 


U(^/  ?"»«■) 

V(%<f 
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(368) 


with 


=  -r»+/3P) 

-  (P-R)  $('*->) 


U)(kv-V 


[£($*) 


+ 


3nA+7k'"-2nx0+Kl) 

tn4  f  -«Y l+kx) 
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(375) 


4/1^  f  K1  -Aa.X(  I  —  u2)  (  I  t  Kl) -t 'Hz  +  ^ 

/Ia^- 


g  is  with  (280)  and  (281) 


184 


+/3P)  ({f; {p) 


-  FQjt, 


m 

be (ff 


It  is 


F(V.  k)  "Z  (  *  )("0 

v»o  v 

£ (<$, x)  ~Z(i)(- 


with  the  recursion  formula 
'  2  V  -  I 


p  2V  -  I  .  v  I  2v-/ 

~T7~  J*v-*(cr)  ~  J7  ^  ^ 

V*'>'r  l 

l  v 

T  f  ~ 


\f*  1 


h  (y)  "  %-p<?)  ~  "jr  **+•?*”¥ 


_  s> . ' 


2>  .  »  „  i  i  .  ,  i,3 

4  2  ¥  4  7  h^if/C^y>  ~  — 


Vt.  3 

fw  -ft  ir 


~  JL  —  JnV'M  >^rocp 
6  4  2  r  ( 


For 


ftL  -  \j^  small  is  Z"  (M )  Kfc)  % 

and 


E(  f  ;  Jt)  ~F(  ?  , 

¥<?-  ^){£(MYm)) 

t  1- 


and  for  null  cp  and  small  k 


£(<f;  fe )  -  F<nk)  ^  “  -  i  k  f 


I  ,  a*  A 
-rr  k  CP 


(388) 


Therefore,  with 


=  -  ,Rf3P)  -  F  mpu 


V ~  I  +  4^  pVy 


•* 

2 


g  can  be  approximated  by 

3  ^ 

»■  ^  ■  T. 

A  2(P-R.) 

therefore,  (368)  with  (389)  and  (390) 


a  _  -ce^p)-^=(-p-g)'fea,4^%£ 

~i  ..  .  i  ^  ~rs — 

1  *  3<r?  7pr^5: 

Setting  X*  -  x,  we  get  the  following  third -order  equation: 

j  (p-g)'^  ,  EzS. 

*  +  *  4  s^Zsf-  +  4 


•2>/z 


pftP)(P-R 

-<fcr-tr. 

Z 


(390) 


(391) 


o  (392) 
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For  small  cp  .the  first  and  third  terms  can  be  ignored,  and  we  get 
max 

X2-^  >  --  (H  -h[3?)  <393> 

as  a  reasonably  good  solution  for  the  eigenvalue  X.  The  eigenvalue 
X  does  not  depend  (or  does  not  strongly  depend)  on  Cp 
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E.  BUCKLING  OF  FINITE  MULT I FIB ERF  IN  A  MATRIX  UNDER  AXIAL  LOAD. 
THREE-DIMENSIONAL  " 

k  FOR  THE  MULTIFIBER  REINFORCED  COMPOSITE  WITH  THE  FINITE 
MATRIX  TREATED  AS  A  THREE-DIMENSIONAL  CYLINDER 


The  displacement  components  of  the  matrix,  in  terms  of  Papkovich 
functions,  in  cylindrical  coordinates,  are  as  follows: 


5  —  ft®**  +  E®**  ~  4{h^)fk 

(395) J 

4  ~ ~  ^  +-?*si«¥) + n]  096) 

where  P  P.  and  P„  are  harmonic  functions  such  that 
o,l,  2 

vij5)  ~  vy* —o  <397) 


y  is  a  Laplacian  operator  in  cylindrical  coordinates;  i.e., 


(398) 
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The  six  stress  components  in  cylindrical  coordinates  are: 

(399) 

C  =i«$5  O'V-  gr,gA'p«*+pS1.«+fi'] 

(401^ 

.i 

^(h?)  (AlT' ^^'4>Pcj 

<£■-*«(«**  ts^S) 

~2V~i?)  (IAS  [A(fieeSH  +  fcj  (403) 

<r;=  gxMH 

'  ^  idfii  sty  + 

\ 
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The  product  solution  form  of  a  harmonic  function  in  cylindrical 
coordinates  is 


(  J*  i  Km  ( Sin  (4*4’)  4  fa  4>)j-  [si*  (o<H 

J  ^  (405)  J 


In  considering  the  boundary  conditions  of  the  present  problem,  we 
choose  the  solutions  of  Papkovich  functions  as  follows: 

f,  ~  [ A, X  («*■*)  +  8,  K, («■*)] 4 5'k c««^) 

^  ^  (406) 


ft  ~  (c,  X,  «*)  AX>,  Kz  K-Vj 

■+  [  ^  lo (*A) tf^  Ko  (% ^)]  5m  ^  (407) 


Then 

P  c«s4>  +  jism< )>  —  [  C(  -E  (<*«  f  J>  K,  («* 

+  F,  KeK4)j  c«54“ -s/V  K^J) 

(409) 

-fs'^  +  =[c(  im 

-F,  KoC<  k)J 

(410) 
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Substitution  of  the  foregoing  equations  (406)  through  (410)  into 
equations  (394)  through  (396)  yields  the  following  solutions  for 
displacements ;  ^ 

|  i  ox,*-)  tx»(  k2  k^) 

tE,  f  Fj  Kc 

J 

-  t  6,  K/faV 


C^Wj+j}  K^A)+E,  j!kx) 


4"  Fj  Fo  j  ^5  4  (411 

40=1?;  fi(A/Ji(^)  +  Bi  K 

-(3-4P1)  [E(IeMfF,K0^A)|- 

•  5m  Sik  (^m^)  (412; 


4)I,Nt6,K,KA) 

trf C,  I liK*)+EIc  KA) 

+  Fi  Ke(«X*4)]  l-CCS<f>  Co(^*) 


The  displacement  components  in  Cartesian  coordinates  are 


=  §*c<rs$  - 

=  ~  s(h^  I /[(A,  I,  fa.*)  t  3,  K,K^ 


freest 

~  %o-p*)  ~k  lA  J-i  ^  $  ^i  ^  A) 


(417) 


The  radial  stress  a??  and  shear  stress  a??  are  found  from  equations 
(399),  (402),  (406) ,  to  407,  ~  1Z 


tF,K/K^] 

+ife)  [-^(A^M+qKMUX) 

-  ( A,.I(  («**)-/-  6,  K^a.)) 

+E,j;(cu)+feK/°u)) 
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CP3  4  s^(p(vC%>) 


(418) 


rr1  =  -S- 

Uyy  "jf  /  — 


+ f  (-zp*-/)  fl 

•  CC3  4s  Si  4  ((*„;&)  (419) 

where 


(420) 
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•'Stub  Sin  (42i) 


-C  [( ' 


•  C^3  ^  CA 3  (V* 


(422) 
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Equations  (414)  through  (422)  contain  six  constants  (A^ ,  B.  ,  ,  E,  , 

and  Fp  which  must  be  determined  by  boundary  conditions.  These  conditions 
include  the  continuity  of  displacements  and  stresses  at  the  interface  and 
will  yield  six  simultaneous  equations  with  six  unknowns  which  are  hardly 
amenable.  For  simplicity,  it  is  then  assumed  that  the  radius  of  the  fiber 
is  much  smaller  than  that  of  the  composite;  the  three  coefficients  for 
modified  Bessel  functions  of  the  second  kind  will  be  automatically  dropped 
because  of  the  finiteness  of  stress  and  displacements  at  the  neighborhood 
of  the  origin. 

The  boundary  conditions  at  interface  r  =  a  to  be  used  are 


5/  -  -  s*  *  t**; 

%  =  ° 

<423> 


Substitution  of  equation  (423)  into  equations  (414),  (416),  and 
(417)  yields 

g(l-p“j  f[a  -T  f  I  A  / 

+ (.2  talWjJc, 

4  -tajfajl  4= r 1 

J  (424) 
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[±  I  K«)  -l!M]At  ife(h^)i(^)-a^K«jjc 
~[aIeW)J  E,  =  o  (425) 

and 

A,  +  C,  3  If^)  t A a_TcK«)— o  (426) 


Solving  these  equations,  we  can  obtain  the  solutions  of  A^,  C^,  and 
Ej^  In  terms  of  5*n> 


The  denominator  of  coefficients  A^,  C^,  and  is  designated  as  A  , 
and  is  written  out  as  below. 


—  2.  j^(/— C^*f 

+  (°^»  a)  +Xf(o(H  4)j  (427: 


The  results  of  the  coefficients  are 

A,  =  ^ 

'  A  ( 

— ■2I*(pt*<ol  (428) 
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c. 


and 


£ltk  i 

A 


«  a)  Jk  K*)IK«)-1p<h 


(429) 


E. 


=  I  -ten*)!/***)  -^-WX)IK«)I^) 


U  +< 


j  jfa,*)  • 


(430) 


The  force  per  unit  length  applied  by  the  matrix  to  the  fiber  due  to 
the  normal  pressure  and  circumferential  shear  at  the  interface  is 

fi~Jo  _ 


lA—cu 


(431) 


By  substitution  of  equations  (419)  and  (421)  into  (431),  and  the 
identities  of  definite  integrals 


JW<H*  =  i  -  ^ 


(432) 


and 


/ 


we  then  have 


i'  —  Sr$r)  (  A,  ( -  <l  i,  M 
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+  £,(<i  J K<C) -H^) 

(433) 


Since  it  was  assumed  that  the  fiber  buckles  on  x-z  plane,  then 


(434) 


Combination  of  equations  (433),  (434),  (414)  and  (425)  through  (427) 
gives  the  unit  interfacial  force  per  unit  lateral  deflection,  k,  due 
to  surface  normal  force  and  circumferential  shear,  and  is  resulted  as 
below: 


O 


s(  K  *>  i?  a) 

+  IS(c/na) 

&  1 

|4  fz(4 


•  I^i***)  I*(<*  *4 


(435) 
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BUCKLING  LOAD  OF  THE  FIBER 

Combination  of  equations  (72)  and  (435)  yields 


r 

I  cx 


IT  E* 


f  [/  4  (j.  -pfy-jjg)  - 

•i;W) 

5~pl)  +3z(/-jfj(2-j^. 

'  I%*a)  I(  K«) 

fK^J  Ik^lW)  j 


where  A  is  defined  in  equation  (427). 
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SMALLEST  BUCKLING  LOAD  AND  BUCKLING  WAVELENGTH  OF  THE  FIBER 

Minimization  of  P  with  respect  to  a  ,  i.e.,  SP  =  0  (437),  gives 
cr  n  cr 

8 <*„,ExXX|“^ c \ -f°4  I(<$t (i -4^)^- 

•iWOxM 

°  '  z 

+ |  (  j- f  j0ka)j,  k«) 

""£-2  C3~4-PE)o<^^  -fo/^3 

•  j  2  ( /-pX>  i3c(„ajUa} 

-  Cp*  (•***')  J?  (<**<) 

-3  a. 

•  [  Ie(°k  *) "  -I,Ko)] 

-  2  [cc,-s>*)  x2M4  (««*) 
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-  (cC**y+  8(l-P) 

•  I,^a)  -“'m 

"f*  A^-Xp  (^H  1 

f 4  [-2  (  4  -J^j  CL  +JoC„P3 

-s  u-jf-x  i-jfpr* 

i  4  j^2  ( 4  -■2-PL~-pL'')  (“<Kp  fK 

+S  i-2-,P^)  ( /-j-^J)  K«T'J- 

•  j^J,  (°<Ha)i'j.Q. fan*) 

-2d«  ' J;^k«) J^WnA)J  l 

r  ^ 

1  ffC  /-^)  fi  «*)I/K  a) 

-[3Jfa4'-i£u-j$(/-j*jJj;%<H4) 

~4{((>  ~^J  K  47-f  s  (i-^c/-o^)]  • 

I,  (**.*) 


f  6  l-p*)  fa*)'']' 

•  Xc  (c**0-)  -Ij  («'««}|* 

aJ^KaJlK*) 

-J-0<H  (Xj  (*k*)J 
-[2(3  cc'-t  m^Ji^)j(2(^) 
-{m  3—q-xO  *«  +°<*<c] -{iffa*) 

•XK^’M 


+  <Ki X  a 

—  o(H  -I,3K<0  oi^l^c^)  X  (<4  5^ 


0 


(438) 


This  equation  determines  the  critical  buckling  wavelength  of  the 
fiber,  -0cr  =  2n /a  i.  After  substituting  the  obtained  value  of 
into  equation  (436),  we  then  get  the  smallest  critical  buckling 
load  of  the  fiber,  (p*r)min 
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APPENDIX  V 


EXPERIMENTAL  INVESTIGATION  OF  SHORT  BORON  FIBER  COMPOSITES 


1 .  MODULUS  TESTS 


The  specimens  consisted  of  a  laminate  comprised  of  8-ply  S-glass 
Scotchply  tape  (50%  by  weight) ,  seven  interspersed  layers  of  approxi¬ 
mately  2-inch-long  boron  fibers  (23%)  by  weight),  and  epoxy  Epon 
828/1031  resin  (27%  by  weight) . 

The  specimens  were  tested  in  compression  and  resulted  in  values 
summarized  in  Table  IIL  For  comparison,  reported  values  of  pure  S-glass 
tape  and  pure  unidirectional  boron  composite  are  listed  in  the  same 
table. 


TABLE  III 

COMPARISON  OF  SHORT  BORON  FIBER  COMPOSITES* 


IN  AXIAL  COMPRESSION 


Load  Application 

Short-Fiber 

Composite 

S -Glass** 

Boron*** 

Max  strength,  psi 

Max  modulus,  10  psi 

133,000 

13.7 

120,000 

8.9 

_ 

220,000 

38.0 

*  Unidirectional 

**  Minnesota  Mining  &  Manufacturing  information 
***  From  Air  Force  sponsored  boron  work  at  Narmco 

In  comparing  the  short-fiber  specimen  with  pure  glass  and  pure 
boron  composites,  it  can  be  said  that  the  values  follow  a  generally 
logical  pattern, which  is  obvious  in  the  graphical  presentation.  Figure  65. 
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2.  FATIGUE  TESTS 

The  fatigue  test  data  of  a  short  boron  fiber  composite  were  com¬ 
pared  with  fatigue  data  of  7075-T6  aluminum  (MIL-HDBK  5)  and  a  composite 
which  was  fabricated  with  a  Scotchply  1002  resin  and  unidirectional 
S-glass  (FPL  Tech.  Rep.  AFML  7R-64-403  October  1964). 

The  properties  utilized  for  these  materials  are  presented  in  the 
following  table. 
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TABLE  V 


MATERIAL  PROPERTIES 


Material 

Density 

#/cu.in. 

Ultimate  Tensile 
Strength,  psi 

Modulus 

10^  psi 

Aluminum  7075-T6 

0.101 

80,000 

S-glass  &  Scotchply 

0.069 

160,000 

6.1 

S-glass  &  short  boron 

0.079 

180,000 

17.2 

Two  panels  were  fabricated  with  unidirectional  short  boron  fibers 
and  unidirectional  S-glass  (single-end  roving) .  The  physical  and 
mechanical  properties  are  given  in  the  following  table. 


TABLE  VI 

PROPERTIES  OF  TESTED  PANELS 


■  ■■  '  *1 

Panel  No. 

Specific 

Boron 

Content 

Vol.% 

Ultimate  Tensile 

Tensile 

Gravity 

Strength,  psi 

Modulus,  psi 

935 

2.11 

11.4 

180,000 

- 

936 

2.26 

17.1 

180,000 

17.2  x  106 

The  fatigue  data  for  each  of  the  above  panels  is  shown  in  the 
following  table. 


TABLE  VII 


FATIGUE  TEST  DATA 


Panel  No. 

Mean  Stress 
%  of  Ultimate 

Alternating  Stress 
%  of  Ultimate 

Cycles  to 
Failure 

935 

30 

±23 

5.2  x  103 

935 

30 

±15 

1.53  x  105 

935 

30 

±10 

6* 

8.65  x  10 

936 

50 

±10 

1.03  x  105 

936 

50 

±8 

2.63  x  105 

936 

50 

±5 

6* 

6.1  x  10 

The  fatigue  data  are  also  presented  in  Figures  66  through  69.  Figure 
66  is  a  modified  Goodman  diagram  for  the  fatigue  life  of  the  short  boron 
fiber  panels,  No.  935  and  No.  936.  The  curves  are  estimates  based  on  the 
data  of  Table  IV  and  previous  experience  with  other  boron  glass  composites. 
Figure  67  presents  S-N  curves  for  the  aluminum  7075-T6,  S-glass  and  Scotch- 
ply,  and  the  short  boron  and  S-glass.  The  curves  are  plotted  as  a  function 
of  the  alternating  stress  (%  of  ultimate)  and  cycles  to  failure  for  various 
mean  stresses.  Figure  68  is  a  comparison  of  the  aluminum,  S-glass  and 
Scotchply,  and  short  boron  and  S-glass  at  0  mean  stress  plotted  as  a  function 
of  actual  alternating  stress  and  number  of  cycles  to  failure.  Figure  69  is 
a  comparison  of  the  aluminum,  S-glass  and  Scotchply,  and  short  boron  and 
S-glass  at  0  mean  stress  plotted  as  a  function  of  number  of  cycles  to 
failure  and  the  specific  alternating  stress.  The  specific  alternating 
stress  is  obtained  by  the  ratio  of  the  alternating  stress  to  the  density 
cf  the  material. 


*  Specimen  did  not  fail. 
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Alternating  Stress,  7.  Ultimate 


Various  Composites  and  Aluminum 


Figure  69.  S-N  Curves  for  Various  Composites  and  Aluminum 
for  0  Mean  Stress  Plotted  as  a  Function  of 
Specific  Alternating  Stress 
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